GEOMETRY AND A PRIORI ESTIMATES FOR FREE BOUNDARY 
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Abstract. In this paper we derive estimates to the free boundary problem for the Euler equation 
with surface tension, and without surface tension provided the Ray leigh- Taylor sign condition holds. 
We prove that as the surface tension tends to zero, when the Ray leigh- Taylor condition is satisfied, 
solutions converge to the Euler flow with zero surface tension. 



1. Introduction 

In this paper we study free boundary problems of the Euler's equation in vacuum: 

(vt + S/ v v = -Vp, i£S] t cR" 



V • v = 0, x G Qf 



where for every i 6 K, v(t, •) is the velocity field of an incompressible inviscid fluid in a moving 
domain (bounded and connected) fit C M n , n > 2, and p(t, •) is the pressure. The boundary of the 
domain Q t moves with the fluid velocity and the pressure at the boundary is given by the surface 
tension, that is 



(BC) 



D 4 = d t + v ■ V is tangent to \J U t C M n+1 
p(t, x) = e 2 K(t, x), x E dtt t , < e < 1 



where n(t,x) is the mean curvature of the boundary d£lt at x £ <9f2j, and Dj is the material 
derivative. This is equivalent to saying the velocity of dttt is given by v ■ N where N is the unit 
normal to d£lt ■ The case e = corresponds to the zero surface tension problem. 

In the presence of surface tension we will derive energy estimates that bound the Sobolev norms 
of the velocity and the boundary. In addition we will show that if the Rayleigh- Taylor sign condition 
is verified, then some of these bounds are independent of e. In this case we conclude that as e — > 
solutions of the problem with surface tension converge to solutions of the zero surface tension 
problem. We do not include the effects of gravity in (E) as it will only contribute lower order terms 
to our estimates. 

The free boundary value problem for (E) has been studied intensively by many authors. In the 
absence of surface tension the earliest mathematical results on the well posedness of the water 
waves problem were given by V. I. Nalimov jNA74j where he considered the irrotational problem 
in 2 dimensions with small data in some Sobolev space (see also H. Yoshihara Y082J. The first 
break through in solving the well posedness for the irrotational problem, no surface tension, for 
general data came in the work of S. J. Wu [WU971 IWU99j who solved the problem in all dimensions. 
For the general problem with no surface tension D. Christodoulou and H. Lindblad |CL00j were 
the first to obtain energy estimates based on the geometry of the moving domain, assuming the 
Rayleigh- Taylor sign condition for rotational flows. H. Lindblad LI05 proved existence of solutions 
for the general problem. In the absence of this condition D. Ebin |EB87j proved that the problem 
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is ill-posed. There is also the work of K. Beyer and M. Giinther PBG98| IBGOOj on well posedness 
which we will comment more on at the end of the introduction. 

For the problem with surface tension H. Yoshihara Y083j, T. Iguchi |IG01| and D. Ambrose 
AM03] solved the well posed ness irrotational problem in 2 dimensions under varying assumptions 
on the initial data. B. Schweizer |SC05j proved existence for the general 3 dimensional problem. 
Recently D. Ambrose and N. Masmoudi AM05a proved that as e — > solutions of the 2 dimensional 
irrotational problem converges to solutions of the zero surface problem by writing the equation in 
terms of the arc length of the fluid boundary. 

There are many other works on this problem we mention the work of D. Lannes LAQ5], T. 
J. Beal , T. Hou, and J. Lowen grub |BHL93| , W. Craig |CR85| . G. Shnider and E. C. Wayne 
|SW02j . M. Ogawa and A. Tani |()T02j . 

During the writing of this manuscript we were informed and received reprints of several related 
work by D. Coutand and S. Shkoller |SC05| . D. Ambrose and N. Mamoudi AM06 , and P. Zhang 
and Z. Zhang [ZZ06 . In D. Coutand and S. Shkoller work they proved local well posedness of the 
general problem using Lagrangian coordinates. 

Our results differs from those mentioned above in that we can obtain the e — > limit for the 
general problem. Our approach to the problem is based on the well known fact that the free 
boundary problem (E, BC) has a Lagrangian formulation given by 

I( u ) = J J ^j-dydt - e 2 J S{u)dt, 

where u(i, •) G T = {<& : Qq — > M n , volume preserving homeomorphisms}, and S(u) is the surface 
area of u(dflt). Critical points of / satisfy 

(E-L) d t u t + q + e 2 S'(u) = 0, 

where q is the Lagrange multiplier due to the constraint u G T. Writing v = ut o u~ l and changing 
to Eulerian coordinates we obtain 

d t u t ^~D t v, q -> Vp V)V = -VA~Hr(DvDv), S'{u)^J = Vk h 

where A -1 is the inverse Laplacian with zero Dirichlet data, k the mean curvature, and Kyi is the 
harmonic extension of k into £lf Thus (E-L) is the Euler equation with the pressure p given by 

P = Pv,v + e 2 K H . 

It is important to note that this derivation splits the pressure into two terms, the first p VyV is the 
Lagrange multiplier, and the second k-h is due to surface tension. Thus these two terms will be 
treated differently in the energy estimates. 

Using this variational derivation, one can interpret the Lagrange multiplier as the second funda- 
mental form of the manifold V C L 2 (£lo,M. n ) and rewrite (E-L) as ' 

%u t + e 2 S\u) = 0. 

where & is the Riemannian connection induced on T by the embedding in L 2 . The above form of 
the equation makes it relatively easy to identify the correct linearized problem 

3>?w + <%{u t ,w)u t + e 2 $i 2 S(u)(w) = 0, w(t, •) G T u(ty) T, 

where M is the curvature tensor of the infinite dimensional manifold T C L 2 . Keeping the highest 
order terms in the above equation we obtain 

(LN) S> 2 w + &o(v)w + e 2 £/ id = lower order terms, 



^Symbols in the Lagrangian description have a bar, e.g. 3), while their Eulerian counterparts do not, 
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where &o{v) is a first order differential operator and stf is a third order differential operator. In 
Eulerian coordinates these terms are given by 

&o(v){w,w)= -V NPv , v \w N\ 2 dS, 

£?(u)(w,w)= f \V T w-N\ 2 dS 

J du(Qo) 

where N is the unit normal and V T is the tangential gradient on the boundary of u(tto). Here once 
again we are led in a natural way to distinguish the two problems in the following manner. 

1) For e > two time derivatives are associated with s$ , which is a positive semi-definite operator 

3 

similar to three spatial differentiation, thus roughly speaking, dt ~ {d x )^ ■ Therefore one may be 
led to believe that the regularity of the Lagrangian coordinates given by dtu = v is | order better 
than v, which reflects the regularizing effect of the surface tension. However this is not true for 
the Lagrangian coordinates since is degenerate, and the regularity improvement of the d£lt is 
geometric and is not reflected in the Lagrangian coordinates system . See Section |2] for examples. 
Thus Eularian coordinates are more suitable to use than Lagrangian coordinates for our estimates. 

2) For e = the leading term involves &o(v) and thus the Rayleigh- Taylor instability may occur 
unless we impose the condition 

(RT) - V N p v , v (t,x) > a > xedn t . 

In this case two time derivatives are associated with 3% which is a positive semi-definite operator 
similar to one spatial differentiation. Thus, dt ~ (d x ) 2 and comments similar to above hold on the 
regularity of dtlf 

3) For e > one can directly obtain nonlinear estimates that depend on e by multiplying (E-L) by 
(® 2 S) k S'. 

4) The control that any power &q(v), with (RT) condition, can give over vector fields is limited by 
the smoothness of the boundary dVtf This fact makes the velocity field v inappropriate vector field 
to estimate because it is smoother than what these operators allow. 

5) Since and 3%q{v) are degenerate for fields which are tangential to the boundary d£lt one needs 
to add the vorticity 00 which controls the rotational part of the velocity which is tangential to the 
boundary. 

These facts imply that a natural energy to control is 
r 1 2 1 

£ = J n 2^ k ~ l ® tJ \ 2 + \\ s/k " 2j \ 2 dx + ^oWK'-V,^- 1 ;) + M^ 3fe _ 1(n) 

where J = V% is less smooth than v and satisfies 

'Si 2 J + &o(v)J + e 2 £/ J = lower order terms, 

In addition to the geometry of T the geometry of dVtt plays a crucial role in the estimates. The 
appearance of k in the surface tension and V k-h in the energy make the study of the geometry of 
the boundary as well as the study of the harmonic extension and Dirichlet-Neumann operators on 
dVlt central to the estimate. In using Lagrangian coordinates these operators may be hidden but 
can not be avoided. 

Based on these estimates one can construct an existence proof using the following iteration 
method. Since the acceleration of the boundary is given by the surface tension plus lower order 
terms, in the first step of the iteration we evolve the boundary using this evolution. In the second 
step of the iteration we establish the evolution of the velocity in the interior. This will appear in a 
forthcoming paper. 

Finally after this work was completed A. Mielke pointed out to the second author that a similar 
geometric approach had been used by K. Beyer and M. Giinther to study the irrotational problem 
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by reducing it to the boundary BG98, BGOO . Indeed they proved local well posedness for star 
shaped domains with surface tension and studied the linearized flow for any irrotational flow. In 
fact they derived the principle part of the curvature of T for irrotational problem which of course 
coincides with our Mq{v) acting on gradient vector fields. 

This paper is organized as follows. In section [21 of the paper we give the intuition behind 
the energy estimates by computing the geometry of V. In section |3] we present the geometric 
computation of the moving boundary. In section |1] we present our energy estimates. In sections [S] 
and we present some examples and basic analytic and geometric calculations. 

Notations. All notations will be defined as they are introduced. In addition a list of symbols will 
be given at the end of the paper for a quick reference. Here we'll present some standard notations 
and conventions used throughout the paper. 

All constants will be denoted by C which is a generic bound depending only on the quantities 
specified in the context. We follow the Einstein convention where we sum upon repeated indices. 

For a domain and x E dflt we denote by N(t, x) the outward unit normal, II the second 
fundamental form where U(w) = V W N E T x dQt for w E T x dQt, and k the mean curvature given 
by the trace of II, i.e., k = trll. The regularity of the domains fit is characterized by the local 
regularity of dfi as graphs. In general, an m-dimensional manifold M C R n is said to be of class C k 
or H s , s > t|, if, locally in linear frames, M can be represented by graphs of C k or H s mappings, 
respectively. For <9fi, throughout this paper we will only use these local graph coordinates in 
orthonormal frames. 

2. The geometry behind the Energy 

In his 1966 seminal paper |Ar66j . V. Arnold pointed out that the Euler equation for an incom- 
pressible inviscid fluid can be viewed as the geodesic equation on the group of volume preserving 
diffeomorphisms. This point of view has been adopted and developed by several authors in their 
work on the Euler equations on fixed domains, such as D. G. Ebin and G. Marsden EM70j, A. 
Shnirelman [Sn85j , and Y. Brenier [Br99 to mention a few. It is this point of view that we adopt 
to explain the motivation for our definition of energy. 

In this section, we heuristically outline our geometric point of view on the free boundary problems 
of the Euler's equation and the intuition leading to the energy estimates in the following two 
sections. Though the discussion in this section are mostly in Lagrangian coordinates, the estimates 
are actually done in Eulerian coordinates in the next two sections. 

2.1. Lagrangian formulation of the problem. One of the fundamental properties of the inviscid 
fluid motion is the law of energy conservation. Multiplying the Euler's equation (E) by v, integrating 
on fit, and using (BC), we obtain the conserved energy Eq: 

(2.1) E = E (Q t ,v(t, ■))=/ ^dx + e 2 [ dS 4 f ^dx + e 2 5(dfi t ). 

Jn t 1 Jdn t Jn t 1 

The main difficulty of these problems is handling the free boundary. A traditional way to avoid 
this difficulty is to consider the Lagrangian coordinates. Let u(t,y), y E fio, be the Lagrangian 
coordinate map solving 

(2.2) ^= v (t,x), x(0) = y, 

then we have v = Ut o u~ l and for any vector field w(t, x), x E fit, it is clear that 

(2.3) D t w = d t w + V v w = (w o u)t o u~ l 
Therefore, the Euler's equation can be rewritten as 

(2.4) u tt = -(Vp) o u, u(0) = idnt, -Ap = tr((£>v) 2 ), p\n t = k, 
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where k is the mean curvature of dQ t . 

Since v(t, •) is divergence free, then u(t, •) is volume preserving. Let 

r = : S7o — » IR™ | ^ is a volume preserving homeomorphism}. 

As a manifold, the tangent space of T is given by divergence free vector fields: 

T$r = {w : Q -» K n | V • «? = 0, where 10 = (u> o $ -1 )}. 

For the remainder of this section we follow the following convention: for any vector field X : 
<J>(Oo) —> hs description in Lagrangian coordinates is given by X = X o With slight abuse of 
notation, we also let = Jq^^q q ) dS, i.e. the surface area of $(Oo). Thus, the energy Eq takes 

the following form in the Lagrangian coordinates: 

(2.5) E = E (u,u t ) = l [ \u t \ 2 dy + e 2 S(u), (u,u t )eTT 

where the volume preserving property of u is used. This conservation of energy suggests: 1) TT 
be endowed with the L 2 metric; and 2) the free boundary problem of the Euler's equation has a 
Lagrangian action 

I(u) = J J ^-dxdt -e 2 j S(u)dt, u{t, •) G T. 

Let @ denote the covariant derivative associated with the metric on T, then a critical path u(t, ■) 
of / satisfies 

(2.6) %u t + e 2 S'(u) = 0. 

In order to verify that the Lagrangian coordinate map u(t, •) of a solution of (E) and (BC) is indeed 
a critical path of /, it is convenient to calculate @ and S' by viewing T as a submanifold of the 
Hilbert space L 2 (n ,M. n ). 

Computing (T^T) ± . For any vector field X : 3>(fi) — > M. n we form the Hodge decomposition 

X = w-Vtp, xp = -A~ 1 V-X, V-w = 0, 

where A -1 is the inverse Laplacian on Q t with zero Dirichlet data. Therefore if $ G T, then 
w = w o <£ g T$T. This implies that normal space of T<j>r at $ is 

(T*r) ± = {-(Vip)o<s>\ip\ d{HQo)) ^o}. 

since the Hodge decomposition is orthogonal in L 2 and is volume preserving. 

Computing Qlf Given a path u(t, ■) G T and v = ut- Suppose w(t, •) G T U ^F, then the covariant 

derivative §>tw and the second fundamental form II u ^(w,v) satisfy 

w t = %w + II u(t) (w,v), fym G T u(t) F, II u(t) (w,v) G (T^T^. 

Let v = Ut o u~ l = v o u~ l and w = w o -u -1 which are in the Eulerian coordinates. Then from the 
Hodge decomposition we have 

(2.7) %w = w t - II(w, v), II(w, v) = -{Vp w , v ) o u, p WyV = -A-Hr(DwDv). 
As we do the estimates in the Eulerian coordinates, sometimes it is more convenient to use 

(2.8) fyw = {£> t w) o u- 1 = T> t w + Vp w , v . 

Computing S'(u). By the variation of surface area formula and for any w G T U T we have 

(2.9) < S'(u),w >l 2 (Q )= / k(w)- 1 dS = / Vnu-wdx 
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where tz is the mean curvature, and k-h is its harmonic extension. Since (Vk^)ou G T u T, we obtain 



(2.10) S'(u) = (Vk h )ou = Jou. 

This vector field J, divergence free on f^, is very important for it connects the free boundary 
Euler's flow with the geometry of dVlt and even of V as we will see later in section 13.21 

Combining 1)2.10)1 and 1)2.7)) with w = ut, we obtain that the equation (|2.6|) for critical paths of 
/ becomes 

(2.11) u u = D 4 w o u = (—Vp V)V — e 2 J) ou, v = u t o u~ x , 

which is equivalent to (|2.4jl . Therefore, the free boundary problem (E) and (BC) is a lagrangian 
system on T given by (j2.fi j) . If e = equation (|2.6|) becomes the geodesic equation on T, which is a 
well-known fact. 

2.2. Linearization. In order to analyze the free boundary problems of the Euler's equation, it is 
natural to start with the linearization. The Lagrangian formulation provides a convenient frame 
work for this purpose. From (|2.6|) . the linearized equation is 

(2.12) ® 2 w + M(u u w)u t + e 2 @ 2 S{u){w) = 0, w(t, •) e T^T, 

where is the curvature tensor of the infinite dimensional manifold T. Below we calculate S) 2 S(u), 
which is viewed as a linear operator on T U T, and M . 

Computing S) 2 S(u). Let g(s, •) be a geodesic on T, g(0) = u. Let w = g s and Q. s = g(s, -)(fio)- 
From 1)2.8)) we have 3> s w = (D s w + Vp w>w ) o g = 0. Differentiating 1)2.9)) . 

& 2 S(u)(w,w) = ^- [ nwNdS. 
ds Jan s 

and substitute the expressions for D S N, D S S, and D s k from 1)3. 1JI . ()3.2j) . and ()3.7)) we obtain 



^^(n)^, to) = / rew- L («ur L + P • w T ) + kD s w ■ N + « • D s iV + ic 1 D > /c 

= / /•cu;~ l (k;u>~ l + I? • w T ) — kV nPw,w — kV w tw ■ N 
Jan s 

+ ur 1 (-AdUsW 1 - w ± \U\ 2 + (P • n)(u; T )) dS, 

where II is the second fundamental form of d£l s . 

Needless to say that this is a very complicated expression for 3) 2 S{u){w,w). We will show that 
Si 2 S{u) is a differential operator and will single out its leading order part. Let us assume that O s 
is a sufficiently smooth domain, then from the trace theorem, 

/a2n/„w- -\ / |V7T„..±|2 



*S{u)(w,w)- I |V w | + kV n Pw,w dS\ < C\w\ H1 ( Qs y 

The third term on the left side is estimated by applying the Divergence Theorem twice, 

/ kV NPw,wdS = I Vh-h ■ ^Pw,w — K-Htv(Dw) 2 dx = / V w w ■ Vn-^dx — / k\7 w w ■ NdS 
Jdn s Jn s Jdn s Jdn s 

Therefore, we obtain 

ffl2c/.,\/.7, „t.\ / |Y7T„.,-L|2 



*S(u)(w,w)- |V w | dS\ < C\w\ L 2 {ms) \w\ H i {ms) < C\w\ Hl{ns) . 
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Much as in the derivation of (|2.9|) . for a general u £T we derive an self-adjoint operator <g/(u) on 
T U T 

*7(u)(w) = (vn(-A du{no) (w\ du(no) )^)) ou 

which satisfies 

J{u){w,w) = [ IV^I 2 dS 

for any w G T U T and w = w o u . In the Eulerian coordinates, takes the form 

^(u)H = VW(-A au(no) (^| au(no) ) ± ), Vu; : u(fi ) -» E n satisfying V • u; = 0. 
Since S> 2 S(u) is self-adjoint, then 

(2.13) ^ 2 S(u) = stf + at most 2nd order diff. operators 



/ 



Computing In the linearized equation (|2.12j) . we need to calculate M{u)(ut,w)ut for a lin- 
earized solution w(t, •). Therefore, we may again assume that u(t, •) is a sufficiently smooth critical 
path of the action I, thus v = v o u is smooth as well, and study the operator M{u){ut, -)ut on w. 

Here we apply a well-known formula in Riemannain geometry formally. For any v, w £ T U T, let 
u = v o u^ 1 and w = w o u -1 

M(u){v,w)v ■ w = II u (v,v) ■ II u (w,w) - II u (v,w) 2 = / Vpv,v^Pw,vi ~ \^Pv,w\ 2 dx. 
For smooth v and w £ L 2 (u(ilo)), clearly \Vp v> v)\l 2 (u(si )) < C\ w \L 2 (u(n t ))- As f° r the term, 
^Pv,v^Pw,w dx = / p VjV tr(Dw) 2 dx = / —V w w-X7p V)V ;dx 

= / (-Vat^)^- 1 ) 2 dS 1 + / D 2 p ViV (w,w) dx. 

Much as in the derivation of (|2.9j) . we derive an self-adjoint operator Mq{v) on T U T, depending on 
u and v, 

M Q (v)(w) = \VH{-V NPv,v{ w \du{n )) L )) ou 

which satisfies 

&o(v)(w,w) = / — VnPv^w^I 2 dS. 
Jdu{n) 

In the Eulerian coordinates, Mq{v) takes the form 

Mq{v){w) = VH(-V 'NPv,vWdu(n )) J ~)> Vw : u ( Q o) -> IK n satisfying V ■ w = 0. 
Therefore, in a very rough sense, 

(2.14) M{u){v,w)v = ^o(^) + bounded operators 

where we used the fact that < &(u)(v, -)v, ■ > is self-adjoint. 



3. The geometry of evolving domains 

Suppose £lt C M n is a family of smooth domains with the parameter t, moving with a smooth 
velocity vector field v(t,x), x € fit. We calculate various quantities related to the evolution of the 
geometry of the domain, which are essential in the energy estimate of the free boundary problem 
of the Euler's equations. 
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3.1. Material derivative D t . For any xq € 0, to , the particle path x(t) is the solution of the ODE: 

x t = v(t,x) x(to) = x 

and the material derivative Dt = dt + V„ is differentiation along the direction of x(t) in the space 
time domain in IR n x E. Clearly, x(t) € dQt if € dQ, to . 

Calculations of D 4 A^ and DtS 1 . At any xo € 9O,t , DtN(to,xo) _L iV(io,a?o) since \N(t,x)\ = 1. 
To derive D i A^(to,^o) ! let r(i) € T x ^dVL t be a solution to the linearized particle path ODE: 

T> t T = V T v r(i ) = r GT xo dQ, to . 

At (t , so), D t iV • r = T> t (N -t)-N • D t r = -(Dv)*(N) • r . Therefore, we have 

(3.1) T>tN= -((Dv)*(N)) T . 
From standard calculations for hypersurfaces 

(3.2) D t dS = (v^~ k + V ■ v T )dS. 

Covariant differentiation D^. For the family of hypersurfaces dftt with the velocity field v we 
define parallel transport along the material line x(t) as follows. Given a tangent vector to € T Xo dClt , 
let r(t) be the solution of the following ODE: 

(3.3) D t r ± T x{t) dn t o D t r = (V rW • N)N , r(t ) = r . 

It is easy to verify that r(i) E T x ^d£lt and that this transport preserves the inner product. 

A natural connection between Td£lt C M n for different £ along the materials lines is provided 
by the above parallel transport which induces the covariant differentiation T)J, the projection 
of Dj in W 1 acting on w £ T x u\dflt C M. n . This covariant differentiation induces the covariant 
differentiations of linear (multilinear) operators on tensor products of TdVLt and T*dCl, which will 
also denoted by T)J . 

Calculation of D^II and D^k. Given r G T Xo d£lt , let r(t) be its parallel transport along the 
material line x(t) which enable us to compute 

(BJU)(t) = -dJ(U(t)) = (D t V T iV) T = (V r D t Ar + V [DtiT] A0 T = P T D t iV + n(([D t , r]) T ) 

From ()3.3|) and (|3.1|) . we have 

(3.4) \D t , t) = D t r - V T (^ + v) = (V T v ■ N)N - V T v = -(V T v) T , 

(3.5) (D 4 T n)(r) = -V T [((Dv)*(N)) T ) - n((V r ^) T ). 

To calculate D^k at (to,^o)> we take an orthonormal frame {ri, . . . ,r n _i} of T Xo dflt and parallel 
transport it into an orthonormal frame along x(t). Thus D^zc = Dt(II(rj) • Tj) = (D t r II)(rj) • 
and ([3,5)1 give slightly different but useful forms for D^k 

(3.6) D t /c = - V ■ ((Dv)*(N)) T - U(n) ■ V n v = -A dn v N-2U- ((D T \ Tdnt )v) 

(3.7) D t « = - A 9n v ± - ^|n| 2 + (V ■ U)(v T ). 

Calculations of commutators involving D<. In the following, we will calculate the commuta- 
tors of Dj with operators fi, AT, and Aqq, to show that they are of lower orders. 

. lr> t ,H]f = A- 1 (2DvD 2 f n + Vf n -Av). 

To start, we write the basic formula for any function f(t,x), x € $7^, 

(3.8) D t V/ = VD t /-(Du)*(V/). 

For the tangential gradient, using V T / = V/ — (Vjv/)^, it is straight forward to obtain 

(3.9) D t T V T / = V T DJ-((IV)*(V T /)) T 
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Let f(t,x), x G dU t , be a smooth function. Recall fn = represents the harmonic extension 
of / into Qf We have 

(3.10) AD t f n = V t Af H + 2Vv ■ D 2 f H + V Av f H = 2Dv ■ D 2 f n + Vfn ■ Aw 
which implies 

D t /« = H(T>tf) + A^ABtfn = H(B t f) + A" 1 ^ • D 2 / w + V/ w • At-). 
Therefore we can write 

(3.11) D t H(f) = H(p t f) + A-\2Dv ■ D 2 f n + Vf H ■ Av). 

[D t , A'^g = A' 1 (2Dv ■ D 2 A~ 1 g + Av ■ VA^g). 
Next, we calculate [D^, A" 1 ]. Let g(t, i), i £ fit be a smooth function and <p = A~ 1 g. From the 
first half of ()3.1f)|) where Af = was not used, 

T) t g = D t A<j) = AD t (j) - 2Dv ■ D 2 <p - Av ■ Vcf>. 

Since T)t4>\dsi t = ^, we obtain 

(3.12) D t A _1 g = A -1 D t + A- 1 (2Z)t> • D 2 A~ 1 g + Aw • VA _1 g) 

• [D tj A/]/ = V N A-\2Dv ■ D 2 f n + Vfn ■ Av) - V f H • Vjvi; - V v t /V • N. 
To calculate the commutator of [D t ,A/], from (|3.1I) . (|3.8|) and (|3.11j) . we have 

D t (V/« • N) = V N T> t f H - Vfn ■ V N v + Vfn • D t iV 

= V 7V [H(D t /) + A^Dt; • L» 2 /w + V/« • Av)} - Vfn • - V v t /W • iV. 

Thus, 

(3.13) D t A/v7) = M{T> t f) + V 7V A- 1 (2Z)t; • D 2 / w + Vfn ■ Av) - Vfn • - V vT/ t; • AT. 

. [Aan t ,D t ]/ = 2P 2 / • ((D T \ TdQt )v) + V T / • A dQt v - K V yTf v ■ N. 

In order to calculate the commutator [Ag^ t ,D t ] at xq G dflf , take an orthonormal frame 
{ri, . . . , r n _i} of T XQ d£lf We first extend this to an orthonormal frame to T x dQ,t for all x G c?fit 
close to xq by parallel transporting {t\, . . . ,r n _i} along geodesies on d£l to starting from x . Par- 
allel transporting them again along the material lines x(t), we obtain an orthonormal frame 
{ri, . . . ,r n _i} of T x dVtt f° r an (ij^O near (io,a?o)- From the standard Riemannian geometry, this 
orthonormal frame satisfies the property that, at (to,xo), T>tj = and fa, 73] = T^nTj — T> Tj Ti = 0, 
which will be used repeatedly. For any smooth function f(t,x) defined on <9fi t , at (to, xq), 

T> t A dn J =T> t (V Tj V Tj f - Vv TjTj f) = V Tj .D t V Tj / + V [T>t>Tj ] V T ,/ - V [Dt)2 , r . T .]/ 

=V TJ V T ,D t / + V Ti V[D ()Tj] / + V [DtiTj] V Tj / - V [Dti p T . ri] / 

=A ant D t / + 2V 2 f{r 3 , \D t , Tj \) + V^. ^^-^v^f- 

For any vector field r(t,x) G T x dQt, it is easy to see that [D t ,r] G T x dQt since (a) r, + u G 
T(Ui5fii) => [D t ,r] G T(U t 5fit) and (b) [D t ,r] = D t r - V r u does not have component. Thus, 
D Tj [Dt,Tj] — [Dt,T> Tj Tj] G TdVtt an d we can drop all the normal components in its calculation. 
Using T> Tj Tj = V Tj Tj + kN and Vtj = at (t ,x ), we obtain at (to,Xo), 

V Tj \D U Tj} - [T> u V Tj Tj\ = {V Tj \Dt,Tj] - D^V^.Tj + kN)) T 



. . ,' ' — 1 # \ ~. — / , s /• \ > — — ! \ ■ ,<: 1 / 1 1 -{- / , 11 / J i ' 1 1 v 1 1 ^ 



=(V Tj D t r i - V Tj V Tj v - VtV^Tj - KD t N) 1 = -(Aan t «) ' + k({Dv)*(N))) 
Therefore, from (|3.4|) . 

(3.14) VtA dQ J = A aQt D t / - 2P 2 / • ((£ t |tm» - V T / • A dQt v + kV^tjV ■ N. 

Calculation of D 2 k. This calculation starts with formula (|3,6|) . Since IT : Tdilt TdQt then 
II • D T \xdn t v = n • V|Tar2i w - Let {ri, . . . ,T n _i} be an orthonormal frame which is the parallel 
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transport of an orthonormal frame T XQ d£lt along the material line x(t) 6 d£lf From (|3.1j) . (j3.5|) . 
(|3.8j) . and (j3.6j) . we have at (io,^o)> 

D 2 k = - B t A aQt v ■ N - A dQt v ■ DtN - 2(D t T (n(ri))) • V n v - 2U(r t ) • D t (V Ti «) 



= - B t A dQt v • AT + A dnt v ■ (Dv)*(N) 1 + 2V n [((Dv)*(N)) ' J • V n v 
(3.15) + 2n((V Tl f) T ) • V n v - 2U(n) • Vd^w - 2n(r 4 ) • V Ti D iV + 2n(r 4 ) • (Dv) 2 {Ti) 

= - T> t A mt v • iV — 211 • (Z^I^D^) + A 9QtV • (Z^)*(iV) T + 2[£> (((Z>v)*(jV))~ 

+ n((D T | Taf ^) T )] • p T |T5n, «) +2n- ((Dv) 2 \ TdQt ) T . 
To compute DjAg^t; • N from (|3,14|) we need the general formula 

V 2 f(r, t') = D 2 f(r, r') - (H(r) • r') Vjv/. 
for any t,t' G T^dO^. Therefore, 



(3.16) 



-B t A dnt v -N = -N- A dQt T> t v + 2iV • D 2 v( n , (V n v) T ) - 2(V N v ■ N)(U ■ (D T \ TaQtQ v)) 
+ N ■ Vv((A dnt v) T ) - K \(Vv)*(N) T \ 2 . 



When v and f2j satisfy the Euler's equation, the expression for D 2 k can be written as 

(3.17) D| k = -N ■ A dnt V t v + 2e 2 n • (D T \ Tdn J) + r 
where we signaled out the important terms in the above equation 

3.2. An important vector field for the water wave problem. Since J = appears in the 
Euler's equation as a force generated by the surface tension and its regularity is closely related to 
that of <9ilj, we will study the temporal evolution of J for the rest of this section. 
Computing D t J and D 2 J. From (|3~d]>. ([3"%]) . (j3~TT|) and the definition of J, 

D t J =VD t K n ~ (Dv)*J = Vn(D t K) + VA" 1 ^ ■ DJ + J ■ Av\ - (Dv)* J 

(3.18) =VH(A aQt v ■ N) - 2VH(U ■ ((D J '| r «i») 

+ VA _1 [2Zh; • DJ + J • A«] - (£>?;)* J. 

From (|3.7j) . a slightly different way to write D^J is 

x W = -VH{A mt v L ) + VH[-^|n| 2 + (V ■ U)(v T )] 

+VA~ 1 [2Dv ■ DJ + J • Av] — (Dv)* J. 

Generally, when the surface tension is of order 0(1), it is sufficient to consider DjJ. However, 
when there is no surface tension or the surface tension converges to 0, we have to calculate D 2 J. 
Differentiating ()3.18j) . we obtain 

D 2 J =D t (VD t % — (Dv)* J) 

(3.20) =S/~D 2 k h - (Dv)*V~D t K H - (Dv)*~D t J - (D t Dv)*J 

=V~D 2 k h - 2(Dv)*~D t J - ((Dv) 2 )*J - (D~D t v - (Dv) 2 )*J. 

Since is harmonic, from (|3.11|) and (|3.18|) . 

B 2 k h =D t [W(D t «) + A^(2Dv -DJ+J- Av)] 

=H(D 2 k) + A~ 1 (2Dv ■ D 2 + Av ■ V)W(D t /c) + T> t A^(2Dv ■ DJ + J ■ Av) 

( ' > " ' ' =H(T> 2 k) + A~ 1 (2Dv ■ D + Av)\D t J - VA' 1 (2Dv ■ D J + J • Av) - (Dv)* J] 

+ T> t A' 1 (2Dv ■ DJ + J-Av) 
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A more explicit expression of Ti.(D 2 n) can be derived from ()3.17|) . The DjA" 1 is another term in 
the above equation that can be explicitly calculated from (j3.12|) and (j3.10j) to write 

B t A- 1 (2Dv ■ DJ + J • Av) 

=A- 1 B t (2Dv ■ DJ + J • Av) + A^^Dv ■ D 2 + Av V)A~ 1 (2Dv -DJ+J- Av)} 

(3.22) =A _1 [2(DD t v - (Dv) 2 ) ■ DJ + 2Dv ■ {D~D t J - DJDv) + T) t J ■ Av 

+ J • (AD t u - VavV - 2i),r J (),,,-) 

+ A- l [{2Dv ■ D 2 + Av ■ V)A- 1 (2Dv -DJ + J- Av)} 

Computing &>tJ and &>ttJ with divergence free v. In the rest of this section, assume V-i> = 0. 
Given any vector field w defined on Qt with V • w = 0, let S>tw denote the divergence free part of 
DiU>. It is easy to calculate that 

(3.23) & t w = B t w + Vp VjW , -Ap v ,w = tr(DvDw), p v ,w\aa t = °- 

As J is divergence free, we will decompose material derivatives of J into the divergence parts 
and gradient parts, i.e. we consider @ t J and 3>tt J, the covariant derivatives defined in ()3.23j) . Then 
we have 

(3.24) ® t J = T) t J + Vp v j. 
For the second order derivative, 

(3.25) ® 2 J = D t @ t J + V PvM j = D 2 J + T>t^Pv,j + Vp v ,mJ- 
where p v> j and p v ^ t j are defined as in (|3.23[) , Using IJM.8JI . (|3.12(l 

D t Vp„,j =VD tftiJ - {Dv)*Vp v>J = VBtA-^Dv ■ DJ) - (Dv)*Vp v ,j 

=VA -1 D t (.Dtj • DJ) + VA- 1 {2Dv ■ D 2 + Av ■ D)p V:J - (Dv)*Vp v ,j 

(3.26) , „ 

=VA' 1 [(D'D t v - (Dv) 2 ) DJ + Dv (D~D t J - DJDv)} 

+ VA- x {2Dv ■ D 2 + Av ■ D)p VtJ - (Dv)*Vp V)J . 

Using the above calculations, we will show in Lemma 14.41 that J satisfies the linearized Euler's 
equation with lower order terms. The estimates on J and J from the linearized Euler's Equation 
will imply the estimates on the geometry of the moving domain and the velocity fields. 

4. Main Results 

In this section, we will derive local energy estimates and prove convergence theorems. We show 
that solutions of (E) with boundary condition (BC) are locally bounded 

(4.1) v(t, •) G H 3k (n t ) and dQ t G H s °, s = 3k or 3k + 1 for e = or > 

where k is an integer satisfying 3A; > t|+1 (equivalently 3k > §+§)• When e > 0, this estimate is ob- 
tained without any additional assumption and it may depend on e. To derive a priori estimates inde- 
pendent of e, we assume the Rayleigh- Taylor sign condition (RT): — V NPv,v(t, x) > a > x G dtlt 
for some constant a. 

Definition of the energies and statements of the theorems. The conserved energy of the 
Euler's equation is given by 

E = E (n,v)= [ hv\ 2 dx + e 2 S{n), 
Jn 2 

where 5(0) = Jq Q dS is the surface area. Higher order energies are based on the linearized Euler 
flow and thus involve the differential operators g/, and Mq(v). 
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Recall that, for any vector field uu 6 H s (£l) with V • w = 0, srf is given by 

(4.2) j^(iy) = -VHA dn w ± . 

srf is a semi-positive definite self-adjoint third order differential operator if uu\gu is not smoother 
than N 6 iI So (Sf2). In fact, it is positive definite acting on the irrotational part W{ r of uu, i.e. 
Wir = VHN^w^, see fTM for details. 

Also recall that M$(v), which depends on Q as well as on a vector field v S H 3k (Q) with V-v = 0, 
is given by 

(4.3) &o(v)(w) = VHH-VnPv^w^ = Vnd-VNA-HriDvDv))™ 1 ). 

By Lemma 16.41 VnPv,v £ H 3k ~ 2 (dQ) and therefore &o(v) is a first order self-adjoint differential 
operator if w is not smoother than Vp VjV . Under the sign assumption (RT), &&o(v) is semi-positive 
definite, and like <e/, it is positive definite on the irrotational part Wi r of uu. 

Let oj v : W 1 — * M. n , often simply written as u for short, represent the curl or vorticity of a vector 
field v defined on Q, i.e. 

u(X) ■ Y = V x v ■ Y - V Y v ■ X 
for any vector X, Y £ W 1 . Viewing uj as a matrix, its entries are ujj = uj(-J^) ■ -J^ = diV^ — djV % . 

Definition 4.1. For any domain Q in H s ° , so = 3k if e = or sq = 3k + 1 if e > 0, and any 

vector field v £ H 3k (Q) with V • v = 0, define the energies E(Cl,v) and Ert(Q,v) , often written 
as E and Ert for short, 

E= ( h^ k - l 9 t J\ 2 + j\^ k 'h\ 2 dx + M|a*-i (n) > E RT = I ^& Q (y)*/ k - 1 J ■ $? k - l Jdx 

where is the divergence free part o/D^ J defined in (|3.23f) . 

&tJ = Di J + Vp v ,j = D t J - VA^triDvDJ) 

Set £ = E + E RT . 

To estimate terms in the energy we need to consider the following type of H s ° neighborhoods of 
fi*, a bounded connected domain in R n , which are bounded in H s for some s > sq. 

Definition 4.2. Let A = A(Q*,sq,s,L,5) be the collection of all domains f2 satisfying 
(Al) there exists a diffeomorphism F : dQ* — > dVt C W 1 , so that \F — idga, \H 3 o(dn,) < 
(A2) i/je mean curvature k of dVL satisfies \n\H s - 2 (dQ.) < ^- 

Fix < 5 < 1 and let L = 1 + 2|k(0,-)| 3fc _s /Q s and A = A{n ,3k - i 3fc - 5,L ,5). By 
Lemma 16.41 and equations (jtj.lUj) and (|6.19|) we have 

(4.4) |^U ( H.(n),ff.-3 ( n)) < C, sG[3,3fc-l] 

(4.5) |^b(v)U(ff'»(n) lJ ff»-i(n)) < C'l^^l^sfc-i,^ < C|u|flafc(n)> s G [l,3fc - 1] 

(4.6) |pu,u|if3fc(Q) < C|(L'u) 2 |^3fe-2('m < C\v\ 2 3k _5 by Sobolev inequalities 

where C is uniform in $7 € Ao- Here used the fact 3k > § + §• The norm H 3k ~s in Q4.6|) is chosen 
for convenience; any norm ff ~ l+a wou i c j wor k with a > \. The next proposition gives bounds 
on the velocity and mean curvature in terms of these energies. 

Proposition 4.1. For Q G Ao with dQ G IP , we /iai>e 

e2 | K llf3fe-i(gQ) < 3E + Coe 2 , I u lfl- 3fe (f7) — Co(-S + -So) 

and, if we also assume (RT), 

\ K \H 3k - 2 (8n) — C*Ert + Co, 
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for some constant C*,Co > 0. C* depends only on a in assumption (RT) and Cq depends only on 
the set Aq. 

The proof of this proposition will be given below. Using this result we will prove the following 
three theorems. The first theorem holds when there is surface tension which makes the regularity 
of dflt better (in H 3k+l ) but the bound on |«|jy3fc-i(gn t ) depends on e, the strength of of the surface 
tension. 

Theorem 4.1. Assume e > and fix 5 > sufficiently small. Then there exists L e > such that, 
if a solution of (E) and (BC) is given by fit with dflt E H 3k+1 and v(t, •) E C{H 3k {fl t )), then there 
exists t* > 0, depending only on \v(0, -)\H 3k (Q t )j L e , and the set Aq, such that, for all t € [0, t*\, 



fit e Ao and |«[^3k-i(an t ) ^ L e, 

E(fl t ,v(t,-)) <2E(fl ,v(0,-)) + C 1 + f P t (E ,E(n t ,,v(t' r )))dt> 

Jo 

where P t () is a polynomial of positive coefficients determined only by e and the set Aq and C e is 
an constant determined only by e, \v(0, Ol^sfc-f/Q y an< ^ ^e se ^ ^o- 

Since the domain is evolving, the above continuity assumption of v in t means that there exists 
an extension of v to [0, T] x W 1 which is continuous in H 3k (W l ). The second theorem holds under 
the assumption (RT) and the estimates are uniform in e £ [0, 1] . As it does not take the advantage 
of the surface tension even if it is present, the bound on the regularity of dflt is only in H 3k . 

Theorem 4.2. Assume e € [0, 1] and (RT) holds. Fix sufficiently small 5 > 0. There exists L > 
such that, if a solution of (E) and (BC) is given by fl t with dfl t G H 3k and v(t, •) G C(H 3k (fl t )), 
then there exists t* > 0, depending only on \v(0,-)\jj3km t ), L, and the set Aq, such that, for all 
t£[0,t*], 

Qt £ Ao and \ K \H 3k - 2 (dn t ) — L, 

(J ' lS) £(n t ,v(t,-))<£(n ,v(o,-))+ f p(£(n t ,,v(t',-)))dt' 



where P(-) is a polynomial of positive coefficients uniform in e, determined by the set Aq. 

An immediate consequence of the above theorem is convergence of solution as the surface tension 
approaches 0. 

Theorem 4.3. Assume (RT) holds. Fix the initial data dV> € H 3k+1 and v{0,-) E H 3k (fL ). 
As e — > 0, subject to a subsequence, the solution of (E) and (BC) with vanishing surface tension 
converges to a solution of (E) and (BC) for e = weakly in the space of dflt G E[ 3k and v(t, •) G H 3k . 

The above convergence of dflt is in the sense of local coordinates and the convergence of v can 
be obtained by using the Lagrangian coordinates u(t,y) which is also in H 3k . We also observe 
that the neighborhood Ao of the domains fl does not have to be centered at flo. Thus, since the 
constants involved in the energy estimates only depend on the neighborhoods and the norm of the 
initial velocity, these estimates provide a basis for a continuation argument local in time. 



Proof of Proposition I4TT1 From the definition of s& and 3?o(v), it is easy to obtain 
e 2 |V T AT(-A 9n Ar) fe - 1 K|| 2(an) = eV fe ~Ml^(n) < 2E, 
I -Vat^IA/X-A^AO^M 2 dS =< ^{v)^ k - l J,^ k - l J >| 2(n) = 2E RT . 

J dft 
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To estimate n in either H 3k ~ 2 (dfl) or H sk ~^(dQ), it is sufficient to use the estimates on functions 

Q7 1 

and operators defined on dQ considered as in H ~2 only. Therefore, the inequality for k in Propo- 
sition ^3 follows from (RT) and the fact that N behaves as a first order derivative (Theorem I6.1[ 
ijfTreji . and JEEl)). 

To estimate v, it is easy to calculate that 

(4.9) Av l = djuj) 

which is part of the energy. Therefore, we only need to show that some boundary data of v is 
controlled by E and the conserved energy Eq. This boundary data of v turns out to be Vat v. Since 
Vjy v = (Dv)*(N) + lo(N) and to is controlled by E, it suffices to estimate v = (Dv)*(N). 

Step 1. Tangential curl loJ of v T . Let u T be the tangential component of v and u>J(x) : T x d£l — > 
T x d£l be defined as 

loI{x)(X) ■ Y = V x v T ■ Y - V Y v 1 ■ X = V x v T ■ Y - V Y v 1 ■ X 

for any x £ d£l and 1,7 £ T x dQ. To obtain a more explicit form of u>J, let X and Y be extended 
to tangent vector fields on a neighborhood of x on d£l by parallel transport along geodesies on <90 
emitting from x. From the definition of u, we have 

luJ(X) ■ Y =V x (Vy v-N)- V X Y ■ v - Vy(V x v N) + VyX ■ v 

=U(X)-V Y v-U(Y)-V x v. 

Therefore, by Sobolev inequalities, there exists C > uniform in f2 G Ao so that 

Icj.Ti , fc 5 <c|noD T w| „, 5 < c\v I , t i 

since 3fc > 5 + |. Again here the norm H 3k ~s is chosen to illustrate that the term is lower order. 
In fact any H 3k ~ a with < a < \ works. 

Step 2. Divergence T> ■ u . At any xq € dQ, let {X\ , . . . , X n _x} be an orthonormal frame of T xo dQ. 
We extend them to orthonormal frames of T x dSl at each a; in a neighborhood of xq in dO, by parallel 
transport along geodesies on d£l emitting from xq. At xq, 

(4.11) V ■ u T = V Xi is T ■ X t = V x .{V Xi v-N) = A dn v-N + (D T \ Tdnt )v ■ U. 

To control the first term on the right side, we use Q3.24JI and (|3.18j) to obtain 

\9tJ + VH(A dn v ■ N)\ H3k - 3{n) < C\v\ H3k _i {Q) 

where we also used Lemma 16.41 for the estimates. To get boundary estimates, note that S>t J 1 N 
and N • VH(Agn v ■ N) are well defined in H 3k ~z (<9$7) even if k = 1 since they are divergence free, 
and thus combining the above inequality with the identity for T> ■ z/ T (|4.11[) . we obtain 

W(V ■ u T ) + (& t J) L \ 3fc r <C\v\ 3k i \V-v\ 2 , h 5 <C(E+\v\ 2 , ll j ). 

5*tep 5. Control of v T . Using the same frame as in Step 2, at xo, 

A dn u T ■ X % =V Xj V x . v T ■ X % = V x , (V x . v T ■ Xi) = V Xj (V Xi u T ■ Xj + ^(Xj) ■ X t ) 
=V Xi (V ■ u T ) + 1l(X u Xj)v T ■ X 3 + (V Xj loJKXj) ■ Xi. 
One can also write 

A 9 n = V T (P • is T ) + Ric((Vt;)*(^) T ) + (V Xj <4){Xj) 
where Ric is the Ricci curvature of dQ. From the estimate on toj and T> ■ v T , 

\v T \\ k a <C(E+\v\\ k ! ) 
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with a uniform constant C > 0. 

Step 4- Normal component u 1 - = Vnv-N of v. This will be estimated by calculating the divergence 
of v in two ways. Recall Nn = (N^, . . . , N^) denotes the harmonic extension of N into CI. Let v 
also be extended to (Dv)*(Nn)- Near any xq E dCl, let {X%, . . . , X n _i} be the orthonormal frame 
of T x dCl constructed above and let X n = Nn- On one hand, at xq, 

V • v =V Xi ((Dv)*(N n )) ■ Xi = V Xl (V Xl v-N n )-N- Dv(V Xi Xi) 
=(V Xl u){Xi) ■ N + Dv ■ DN n 

where the first term in the last line follows from (|4.9jl . Therefore 

|V - v - (V Xi u)( Xl ) ■ N\ H3k _ km) < C\v\ H3k _ h (fl) . 

On the other hand, by decomposing v into the tangential and normal parts, one may calculate V • v 
alternatively 

V • v = V ■ u T + ku 1 - + V N u ■ N = T> -u T + kv l + V N (V Nn v ■ N H ) - N ■ Dv(M(N)) 
which along with the previous identity implies 

\V ■ v T + V N {V Nn v-N H )- {V Xi u){Xi) ■ N\ H3k _& {m) < C\v\ H3k _i (n) . 

Since u> is controlled by E and T> ■ v has been estimated in Step 2, we have 

|V«(V J v„».JV„)|^. }(8s!) <C(f ; + H^. J([i) ) 

The following decomposition trick on dCl has been used many times in the basic estimates in 
Sections El and |31 

Vjv(Vjv h v ■ N H ) =M{V N v N) + Vj\rA _1 (Vj\r„Au • N n + 2tr((DN H )*DvDN H ) 
+ 2D 2 v(N n ) ■ DN H + 2D M^, ^) • N H ). 
Using (|4,9|) again, we have with a uniform C > 0, 

From Step 3 and Step 4 above , we have 

M 2 3fc 3 <C(E+\v\\ k ! ) 

which implies the same estimate of the boundary data V x v on dCl. Combining it with the Poisson 
equation (|4.9|) . we obtain 

Mh m <C(E + \v\l 3k _ hQ) ). 
Since \v\ 3fc _i < /3|u|jj3fc/n\ + Cr\v\i2(q\, proposition B~T1 follows immediately. □ 

In the proof of Theorem 14. II and Theorem 14.21 we will need the following lemmas. 
Lemma 4.2. For any CI E Ao with k E H s (dCl), s E [3k — |, 3A; — |] ; we /iai>e 

|n|^(en t ) + \N\ H *+i(dn t ) - ^(1 + M-ff s (an t )) 
/or some C > uniform in Cl E Aq. 
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Proof. We only need to prove the estimates for IT. This is obvious if n = 2. For n > 3, we use 
identity (fH3|) : 

-A an n = -v 2 k + (|n| 2 / - kU)u. 

Since f2 G Ao, we have 

l n '// 3fc -§(9c t ) + ^H Sk 'i(an t ) ~ C 

and it implies 

2 ' 2 



n| 2 / — «n|#»i (gQ t ) < C, where s± = min{6/c — 5 — ,3k — -} 



as long as 3k - § ^ 2=1. Therefore, for s 2 G [3A; - |,3fc - |]\{2fl}, 

|(|n| 2 /- KlI)II| HS 3 (ant) < C\U\ H s 2{mt) s 3 =min{si,si + s 2 - ^y^ 1 }- 

Since S3 + 2 — s 2 > min{l, 3fc — 1 — ^} > and si + 2 > 3k — |, the estimate of II can be improved to 
H s (d£lt) by bootstrap on Ag^. The exceptional cases of the indices can be handled similarly. □ 

Although this regularity of II in terms of k, follows directly from Proposition 16. 3( the point of 
this lemma is that the constant C > depends only on Ao, i.e., f2 G H 3k ~z . This is also the point 
of the following lemma. 

Lemma 4.3. For fl G A with dtl G H 3k and v G H 3k (Vt), we have 
for some C > uniform in f2 G Ao- 

Proof. The idea of the proof is to use p v , v \dn = and the identity / = 7i(f\dn) + A _1 A/ for any 
/ : 0, — > M. On the one hand, notice on <9f2, 

VjvVjv w p„,„ = Af(Af) • Vp U)U + D 2 p VtV (N, N) = M(N) • Vp V;V + Ap v>v - kV nPv>v , 

which, along Lemma 14.21 and ()4.6|) . implies 

\^N^N n Pv,v\H 3k - 2 (dCl) < C(l + \K\ H 3k-2( d ty)\v\ 2 H3k ( n y 

On the other hand, using Lemma 14.21 and l|4.6j) as well, 

\^N n Pv,v\ H s k ^ (n) =1 ™ • Vtr(D V ) 2 + 2D 2 p^ • £)jV w | H3fc _5 

<C(1 + |K|^3fc-2( 9 Q))|w| 2 ^ 3 fc( n) - 

Therefore, we obtain the estimates on V n h Pv,v 

The estimate of D 2 p v ^ v is also achieved similarly. Firstly, 

|ADV„| H 3 fc _r (n) = iZ^trpt,) 2 !^^ < C| W | H3fe(n) . 

For the boundary value of D 2 p VjV , we first consider D 2 p VjV (X, X) at x G 80, with A G T x dQ.. As 
usual, extend J to a vector fields in a neighborhood of x on cTi by parallel transporting X along 
geodesies emitting from x. Thus, 

D 2 Pv,v(X, X) = VxVxPv,v ~ ^v x xPv,v = V N p VjV U(X,X). 

Also, we have 

D 2 Pv,v(N,N) = Ap v , v - kVnPv,v = -ti(Dv) 2 - KV N p V)V 

D 2 p VjV (N,X) = Vx^NPv,v ~ ^V x NPv,v = Vx^NPv,v 

Therefore, from the above estimate on V]\r n p vv , we obtain the estimate on D 2 p vv . □ 



WATER WAVES 17 

The following lemma is the most important observation of this paper, which states that J is a 
solution of the linearized Euler's equation up to lower order terms. For the rest of this section, let Q 
denote a generic positive polynomial in \v\ H 3k(QA, \^\H 3k - 2 (dQ t )' anc ^ e \ K \H 3k - 1 (dfi t ) with coefficients 
depending only on the set A . 

Lemma 4.4. Suppose a solution of the Euler's equation is given by S7f £ Ao with d£lt £ H s ° , 
sq = 3k if e = or so = 3k + 1 if e > 0, and v(t, •) € H 3k (Q t ), then we have 

\3> 2 t J + &o(v)J + e 2 £/J\ H3k - 3{nt) < (1 + \@ t J\ H3k - HQt) )Q. 

Proof. We start our proof by two simple observations. First we recall that for any vector field w 
defined on fij satisfying V • w = 0, we write p VjW = 

(4.12) \Pv,w\H°+i(n t ) < c \ v \HM(n t )\ w \H°(n t )i s G l°> 3k ~ ^ 

where C > is uniform in f2 6 Ao- Second the first order derivative $>tJ is lower order since 
by (jg]5H) we have 

(4.13) \@ t J--D t J\ H s (Qt) <Q s€[0,3fc-l], 
and by (|3~TS1) and (|3~5f) 

(4-14) l D ^l^-3(,9f7 t ) < CklifSfc^). 

To verify the lemma we consider the expression for S> 2 J given in H3.25j) and keep the least regular 
terms. Thus 

\& 2 J - T) 2 J|i/3fe-3( nt ) < Q + |DtVp„,j|.ff3fc-3(n t ) 
is given by (j3.26ft and can be estimated using (|4.12|) . lemma 14.31 and Euler's equation 
l D tV^,j| H 3fc-3 (nt) <Q + C\DB t v ■ DJ\ H 3k-4 {Qt ) < Q. 
To estimate T) 2 J given in (|3.2()j) . we use (j4,fi[) and Euler's equation to obtain 

|D?J- VT> 2 t K n \ H3k -3 m <Q. 
To estimate VD 2 k h we use (|3~2Tj) and (j3~2^|) 

\VB 2 k h - VH(D 2 K)\ H3k - 3{nt ) < Q- 
Combine these inequalities, we obtain 

\& 2 J-Vn(B 2 K)\ H3k - 3(nt) <Q. 
The term T) 2 k has been calculated explicitly in (|3,17|) which yields, 

(4.15) \3> 2 J + S/H(A dnt B t v ■ N - 2e 2 U ■ (D T \ Tdnt J))\ H3k - 3{Qt) < Q. 
To deal with Ag^D^t; • N we use Euler's equation to obtain that on the boundary 

A dQt B t v ■ N - 2e 2 n • (D T \ TXh J)) = -N ■ A d n t (V Pv , v ) - e 2 A dn M( K ) + e 2 J ■ A dUt N. 
The last term can be bounded by the identity AgQ t N = —\T1\ 2 N + V t k, and Lemma 14.21 
\^-^n t N\ H3k _ kmt) = \e 2 \V^\ 2 -e 2 M( K )\U\ 2 \ H3k _ kmt) < Q 

Substituting the above into (|4.15|) and Lemma 14.21 

(4.16) \@ 2 J- VH(N ■ A dnt (y PvjV )) + e 2 afJ\ H3k - 3{Ut) < Q. 
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As we are very close to the final desired form, the second term on the above left side has to be 
related to Using formula (|6.13j) . on dQt, we have 

- N ■ A dQt Vp V)V = -N • AVp V)V + kN ■ V N Vp V)V + N • D 2 {Vp v , v )(N, N) 



=N ■ V{tr(Dv) 2 ) + V N (k h Vn h Pv,v + D 2 p V)V (N n , N n )) - VnPv.vJ- 1 

- nVp %v ■ V N N n - 2D 2 p V)V (N, \7 N N n ) 

Keeping the least regular terms and using lemma 1131 implies 

\- N ■ A d u t Vpv,v + Vatp^J" 1 ! 3fc _s < Q 
(4.17) ( } 

+ |Vat {K H V Nn p v , v + D 2 p VtV (N n ,N n )) \ H 3k-^ mt) 

Let / = K H V Nn p ViV + D 2 p v ^ v (N n , N n ) defined on O t , since p v>v \au t = then 

\f\an t \ H 3 k -3 (mt) = \&Pv,v ~ A dQtPv>v \ H3h _s idnt) = MDv) 2 \ H3k ^ {dnt) < C\v\ 2 H3k{nt) . 
Moreover it is easy to check from Lemma 14.31 

\Af\ H 3k-3(Q t J < Q. 

which implies |Vjv/| 3k- 5 ^ Q- Therefore (|4.17|) together with the definition of and the 

half derivative behavior of VTL implies 

\&oJ + VW(iV • A dnt (V Pv>v ))\ \ H 3k-3 iQt) < Q 

which together with (|4.16f) concludes the estimate in the statement of the lemma. □ 

Proof of Theorems 14.11 and 14.21 To prove Theorem 14.11 in addition to Proposition 14.11 
we need the following: a) the estimates on the Lagrangian coordinates map and consequently 

K G H 3k ~2(dflt), b) estimates on uj = Dv — (Dv)* , and c) commutators involving D^, mostly have 
been done in Section |31 In the following all constant C > will be determined only by the set Ao- 

Estimate of the Lagrangian coordinate map u(t,y). From our assumption on v, the ODE ut(t,y) = 
v(t,u(t,y)) solving u is well-posed. Since u(t, ■) : S7o — ^ fJ* is volume preserving and 3k > ^ + 1, it 
is easy to derive, for any s G [0, 3k], and / G H s (Q t ) 

|/ ou(V)Lff*(fi ) ^ C \f\H^n t )Ht,-)\ s H3k{ Q y 

where C > depends only on s. The proof follows simply from induction and interpolation. By 
duality, for s G [0, 3k], 

\f ° u(t, -)|ff-(no) < C\f\ H s (nt) \u(t, •r 1 !^^). 

Therefore, 

(4-18) \u(t, •) - I\ H s H a o) <C I \v{t', ■)\ H 3 Hnt) \u(t' , -)\% k[no) dt', 

J 

where C > depends only on n and k. Let /i > be a positive large number to be specified later, 

(4.19) t = sup{t | \v(t',-)\ H3k[nt) < fi, W G [0,t]}, 

We have t > due to the continuity of v (t, •) in H 3k (fl t )- Then, for all t G [0,t ], 

\u(t, •) - I\ H 3k {Qo) < \i f \u(t', -)|^3fc (no) dt'. 

Therefore, from ODE estimates, there exists t\ > and C2 > which depend only on \v(0, -)\H 3k (n t ) 
such that, for all < t < min{io,ii}, 

(4.20) Ht,-)-I\ H 3 km <C 2 t. 
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It implies the mean curvature estimate, for all < t < min{io,£i}, 

(4.21) k(*,OI ™ s <l«(0,-)l ->> s + C 3 i. 

Here it is easy to see from local coordinates constructed in Section H3 that C3 is determined only 
by \v(0, -)lH 3fc (n t ) an d the set Ao- Therefore, there exists ti > determined only by |v(0, -)l.ff 3fc (fi t ) 
and the set Ao such that Qj € Ao for < t < min{io,i2}- 

Evolution of the curl uj = Dv — (Dv)* . From equations (E) and (|3.8[) . we have 

D t w = DDjii - (DD t u)* + ((£>v)*) 2 - (£>v) 2 = ((Dv)*) 2 - (Dv) 2 = -(Dv)*lv - loDv. 

It is clear how to obtain the estimate of uj in terms of v. differentiating the above equation 3A; — 1 
times, multiplying it by Z) 3fc-1 iij and integrating it on f^, we have 

(4.22) — / |w|^ 3fc _ 1(n4) ds < C|u| H 3fc ( n t )M^3fc-i (nt ). 

j \if 

The commutator involving T>f First, from (|3.14|) and (|3.13() . it is easy to verify that, for any 
function / defined on dflt, 

7 1 

(4.23) \[Vt,Adn t ]\L(H s i(dn t ),H ! >i- 2 (dnt)) - c \ v \H^ k (n t ) s i G -3fe,3A;- -], 



[ D i,A/]|i(^2(gn t ),i^2-i(9n t )) < C\v\ H 3k(n t ) s 2 G - -]. 



■2 2 J 

si 

To extend the range of S2, we use the weak form of [Dt,A/]: 

/ g[D t ,N]fdS= j g((Dv)-(DvY)(V T f)-N + gN(f)V N vNdS 
Jdn t Jdn t 

+ / gnVfn ■ Av - 2Dv(Vf H ) ■ Vgn + Vg H • VA D \2Dv ■ D 2 f n + V/« • A«)dx. 

To conclude that the above estimate for [Dt,7V] holds for S2 = \- By interpolation, 

( 4 - 24 ) ltDt»^]U(ff«(ait),ff-i(8n«)) < CMff 3 *(n t )> s G [^,3fc - -]. 

Evolution of E: first look. Recall the expression of -Ert and E written as E = I\ + /2 : 



.2 , ,.2 



E RT = \< (@o(v))^ k ~ l J^ k - l J >mn t) = \ [ -VNp^K-MAenJ^JV^l 2 dS. 



2 

Also recall that 



—dS=(D-v T + Kv ± )dS. 



Since 

V • v\ d Q t = V-v T + kv ± + V N v ■ N = 

then 

Ikv 1 - + V ■ v T \ , fc 3 = IVatu • JVI tt 3, < CIvlrTSfcro^ 
and thus ^dS would not complicate the estimates since 3k > § + §• 



I: H/i-e 2 <^J,^-^ ( J> i2(S]() | <Q. 
To prove the inequality I, we use ()4.23|) and (|4.24|) to obtain 

\j t h - e 2 < (—Agn t Af) 2k ~ 1 K,Af~DtK >man t ) I < Q, 
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and from Q4.12JI and (|4.13|) . we have 

(4.25) \%J - VH(B t K)\ H3k _ knt) < C\v\ H s Hnt) \J\ HSk _ hnty 

It implies the estimate I for 4rli- 

II: p 2 - < ^ k ~ l $ t J^ k ~ l ® 2 J >I?m | < Q. 
If k = 1, which may happen when n = 2, 3, 

j t h =< @ t J,@?J >V(a t ), 
where we used the fact < @ t J, {& 2 - T) t @ t )J > L2 ^= 0. If k > 1, 

h = \ I • (-A dnt )(-MA dnt ) 2k - 3 (@ t J) ± dS. 

From (|4~2"3"|) and (|03J), we obtain 

4/ 2 - I T> t {{® t J) L ) ■ (-A^ t )(-A^A ant ) 2fc - 3 (^J) ± dS\ < Q. 
at Jdn t 

On on t , 

T> t ($ t J ■ N) = (Dt@tJ) ■ N - V {mJ) TV ■ N = (® 2 J) ■ N + V nM j - V { o JtJ) TV ■ N, 
which implies, along with (|4.12|) . the estimate II for I 2 . 

Ill: \iE RT - < ^-^o^jy^^fJ >mn t ) I < Q- 

In general, for any function f(t, •) defined on dQt with Q t € Ao satisfying J Qn fdS = 0, we have 

^ / -V^,,,/ 2 ^ = [ -VNp^finv 1 - + V-v T )- 2V NPv , v fr> t f - f 2 -D t (V NPv , v )dS 
at Jdn t Jdn t 

Therefore, we obtain from (|4,6|) . 

4 / -^NPv,vfdS - [ -2V NPv , v fB t f - f 2 ~D t (V NPv>v ) dS\ < Q. 

at JdQt JdU t 
Commuting D t with V and A -1 by and (|3~T2|) . 

D t (JV • V Pv , v ) = N ■ {VT> t Pv,v ~ (Dv)*Vp VtV ) 
= ~ VnPv,vVnv ■ N + V N &- l (2Dv ■ D 2 PviV + Av ■ V Pv , v ) - V N A- l (T> t tr{Dv) 2 ), 



and using Euler's equation (|2.11|) to get 



^D t tv(Dv) 2 = -tv(Dvf - D 2 p VjV ■ Dv - e 2 D 2 n H ■ Dv. 



Therefore, 



4 I -V NPv , v f 2 dS - 2 [ -V NPv , v fB t f - e 2 f 2 V N A- 1 (D 2 K H ■ Dv)dS\ < Q 

at JdQt JdQt 



Substituting / = (-J\FA d u t ) k - 1 J\f(K), we obtain 
d 

■Ert — 

ion 



Ert - I -V NPv ^(-NA a ^ t ) k - l M{K) ■ "D t {-NA ailt f' l M{K) dS\ < Q. 
JdQt 



From (14^1) . (l4"24l . and 

, d 



j t E RT - < ®z(v)^ k - l J^ k - l ® t J > L * m | < Q. 
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In order to apply lemma |4~11 to the estimate III, we need to estimate 

I 3 4 < /-i^jy^-^^j >L2{Qt) _ < ^{v)^~ x J^ k - x %J > i2(Ct) 

{@ t J) x • (-A Sf2t ^) 2fe - 2 (-V^ lV ^) 



an* 

- (^J) x • (-A ant AO fe_1 ((-Vjv^M-Aent-V)*"^) d& 

Our strategy will be to move and the multiplication operator by —VnPv,v in the first integrand 
by commuting them with M and Agn t . Thus, we need to estimate [AoQ t , V nPv,v] ; [Af, Vjyjfyu] , 
and [Ag^tjA/]. Using Lemma l4~3| (|6.12j) . and (|6.16j) . we have, 

Tl — 1 

|[Aao t , VNPv,v}\ LiHS{mt)tH s-§ {dUt)) <Q s G (3 - 3k, — ) 

1 n — 1 

I [AT, Vatp„J| i <Q s G [-. ) 

|[Aan t ,A/]| , s . s, <C s G (- - 3k, 3k - 1). 

Therefore we obtain I/3I < Q, which implies the estimate III on 4tErt. 
Evolution of 8. Combining the estimates I, II, III, and Lemma 14.41 we obtain 

(4.26) \^£\<Q, where Q = Q(\v\ H3k{nt) , |K| H s*-2 (ait )> e|«| H 3 fe -i (9nt) ) 

is a polynomial with positive coefficients that depend only on Ao- This inequality on [0, min{to, £2}] 
where to is defined in (|4.19j) and £2 is determined only by \v(0, ')\jj^(iit) and the set Ao. 

Proof of Theorem \4.^\ Assume (RT) holds. From Proposition 14. 11 (|4.26j) implies (|4,8[) . In addition, 
by choosing [i large enough compared to the initial data, to is bounded below by a constant t* > 
depending only on |<y(0, -)I.H" 3& (fio) an d the set Aq. Theorem 14.21 follows immediately. 



Proof of Theorem \4-l\ Assume e > 0. From inequality (|4,26|) and Proposition 14. 1( we obtain 

(4.27) E(t) - E(0) + E RT (t) - E RT (0) < J QM^*^), l"lfl»*-i(«V)) dt' , 

where we use Q e to represent the dependence Q on e. From proposition 16.61 and (|4.6[) we have 

\E RT \ < C|V A rp i , il ,| i oo (9Ct) |K|^3fc-2 (9Qt) < Cbl^ 3 fc-| (nt) l K l//3fc-2 (9Ct) . 

Interpolating v between H 3k (Q t ) and H 3k ~i(fl t ) and k between H 3k ~ 1 (dQt) and H 3k ~i (d£lt), we 
obtain from Proposition 14.11 

for some integer m > where the constant C\ , which include |/e| 3fe _ 5 , is determined only by Eq 

and the set Ao. Since D t u = — Vp VjV — e 2 J is controlled by E in H 3k ~2 (f2 t ) due to Proposition 14. 11 
we can use the Lagrangian coordinate map u(i, •) to estimate \v(t, -)| ff3fe _3 ^ — |v(0, Ol^fc-f / n \- 
Through a similar procedure of the derivation of (j4.20j) and using Proposition I4.1[ there exists 
£3 > 0, depending only on \v(Q, -)\H 3k (n t ) an d the set Ao so that for < t < min{to,t3}, 

-t 



'^•)'^| ( , t) -|^-)|- fc -3 (no) 



< / Qedt' 

Jo 
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for some polynomial Q t with positive coefficients. Therefore, 

tfjff <itf + Ci(l + |«(0,.)|™ )+ f Qedt' <]-E + C x + /*Q e < 

where Ci is determined only by |u(0, -)l and the set Ao. Thus 

S(nt,t;(t s .))<2^(no,«(0 s -)) + C e + / Qerf*'. 

Jo 

By inserting the above inequality into (|4.27j) and using Droposition l4.il we obtain Q4.7|) . By choosing 
large enough compared to the initial data, Theorem 14. II follows. □ 

5. Examples of Lagrangian coordinate maps less smooth than 8Q t 

In Section|3J we established a priori estimates of the free boundary Euler's equation. In particular, 
the estimates indicate that 8£lt is ^, if e = 0, or |, if e > 0, derivative smoother than the velocity 
fields v\g£i t . This is an improvement compared with the regularity directly given by the ODE 
defining the Lagrangian coordinate maps. It is very natural to guess that the Lagrangian coordinate 
maps might be smoother as well. However, the following examples show that the Lagrangian 
coordinate maps is only as smooth as the velocity fields. 

5.1. Case 1: with surface tension. This is a relatively easy case for the construction of the 
example since we need not to worry about the sign assumption (RT): —VnPv,v > 0. The example 
is given for n = 2 and Qt = ^(l), the 2-dimensional open unit ball. In the polar coordinate, it is 
easy to verify that 



8 

v{t,r,9) = e(r) — , p(t,r,9) 



J r'e(r') 2 dr', supp(O) CC B a (l)\{0} 



is a stationary solution of (E) and (BC) for e = 1. The Lagrangian coordinate map 

u(t,ro,9 ) = {r ,6o + t9{r )) 

is only as smooth as v. 

5.2. Case 2: without surface tension. We will construct an example in R 2 again, which satisfies 
the sign condition —VnP > 0. Consider the domain and the vector field in the form 

fi t = {(r,0)|ri(*)<r<^(t)}, r[ = ^-, r' 2 = ^ 

r\ r 2 

, nS aUt) 8 ,8 
with the functions a\(t) and ®{t,r) to be determined. In the polar coordinate 

t-r ^ r, *-r 9 _ 8 18 _ 8 8 

dr Or ar 80 oe Or r 80 oe 00 or 

One may calculate 

Dv( ® ) — ® + (O + ®) ® Dv( ® ) — rO ^ + ® 

dr r 2 dr r r 80' 80 dr r 2 89 

Thus, it is clear that the above given form ensures v is divergence free: 

,8,8 ,18,18 
V-v = Dv( ]fr ). irr+ D v{ --).-- = 0. 
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Moreover, the above calculation implies 

_ „ ,a'At) ai{t) 2 d _ a>i(t) „ 2ai(i)^, d 

-Ap =tr((Vv) 2 ) = - 2r60 r - 2G 2 = - -(r 2 e 2 ) r . 

Therefore, p must be radially symmetric, i.e. p = p(t,r). It is straight forward to calculate 

„ d ,ai(t) 2 ~ a 2 (t) N d 

or r or 

for some function a^{t). Prom the boundary condition p(t, r±(t)) = p(t,r2(t)), we need for the 
Euler's equation, 



(5 1) M = = ^ r fX x K^k ^ rQ2dr 

\ Dt e + = e, + ^e r + = o. 

Let the Lagrangian coordinate map be r = r(t,r ,9 ) and 9 = 9(t,r ,9o), with r G [^(O), r 2 (0)]. 
Due to the symmetry of the vector field, we have r = r(t, To) and 9 = 9q + 9i(t, tq) which satisfy 

d t r = ^ ^! = e(t,r). 
r 

It is easy to see 

r 2 =r 2 + 2,4(t), A(t)= I a x {t')dt' . 

Jo 

Thus, the system (|5.1(l is equivalent to 



d tt ei(t,r ) = --M^i = -(^log(r 2 + 2A))d t 9 l . 



(5.2) 

1 t ),,!!, I I IV, I — — 

The system (|5.2j) can be viewed as an ODE system for (A, #i(-)) on an open set in the Banach space 
I = ExC°([ri(0),r 2 (0)]). Therefore the unique existence of solution to this system is guaranteed. 
One can solve for 

ft r 2 r 2 _ n a 

9(t,r ) = 9 + Q(0,r ) dt' 0(t,r) = —9(0, vV 2 - 24). 

JO r o + *si[T) r 

Thus, it is clear that the Lagrangian coordinate map is not smoother than the velocity field. Finally, 
one notices that small O would ensure the sign condition p r (t,r2(t)) < and p r (t,ri(t)) > since 
small O would make v = -^--§p + ©(^ r )^ a small perturbation of the irrotational solution 
v = ^p^JL, with a slightly different Si, which satisfies the sign condition. 



6. Appendix I: Basic estimates 

In free boundary problems, it often happens that the moving domain is of class H s and 
moves with an H s ° velocity field with sq < s. Moreover, the estimates usually involve functions 
and vector fields defined on Oj and d£lt- Therefore, in this section, we consider collections A of 
domains f2 which are H so close to some reference domain and bounded in the H s class in some 
sense to be defined rigorously. We will outline some basic estimates on functions defined on and 
dQ and some related operators. Through tedious derivation, these estimates will be guaranteed to 
be uniform for all 0, G A. 
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6.1. Sobolev norms. Let C M n be a bounded connected domain, viewing H s (£l), s > 0, as a 
quotient space of H s (W n ), define the norm 

|<7|^ ( n)=inf{|GV (R n) : G G H s (M. n ), G\q = g} 

where | • is defined through the Fourier transform. As usual, for s > 0, Hq(£1) represents 

the closure of C^°(0) in H S (Q) and H~ S (Q) is isometric to (Hq(Q))*. It is important to note 
that with this definition of H s norm the constants in Sobolev embedding (H s — > L p or C a ) are 
independent of Q. The relationship between this definition of H s norm and the standard definition 
will be explored later on page 1281 

C 1 n H 2 boundary dCl. To consider functions defined on <9f2, let £1 C R n be a bounded connected 
domain with <9f2 of class H 2 C^C 1 . Consider the local graph coordinates of dQ, in orthonormal frames. 
When two coordinate charts of this type overlap, it is easy to verify that the transition map between 
these two local coordinate maps is also of C 1 n H 2 . Therefore, on dCl, the definitions of spaces 
C 1 (dQ) PiH 2 (dft) of scalar functions and C° (d^l) D H 1 (d^l) of (fc,Z)-type tensors, though defined in 
local coordinates, are independent of the choice of local coordinates. The Christofell symbols and 
the usual geometric quantities of the hypersurface <9f2, such as the second fundamental form and 
mean curvature are well-defined in L 2 (dil) and the sectional curvature is in L 1 (9il), for it is like 
the square of the second fundamental form. As these will be referred to later, we give the explicit 
formula in local coordinates here. Let {ei, . . . , e n } be an orthonormal frame and (x 1 , . . . , x n ) be 
the coordinates associated with this frame. Suppose Q locally is given by x n > fix 1 , 
with / £ H 2 n C 1 , then using (x 1 , . . . , x n 1 ) as the local coordinates, we have 

The Christofell symbols T% = MSoL ; 

3 s/T+WW 

d d dif 
The second fundamental form n(— —) • — — ? — 



, x n 



'dx iJ dxi y/l + |V/| 2 
(6.1) Mean curvature « = -d, I - djf ] = ^ r + < ): fi) 1 



yi+iv/iv (i+iv/i 2 )5 (i+iv/p) 

d d d d N _ dufdjjf - (dijf) 2 



3 ' 
2 



Sectional curvature TZ{— V ^, _) _ i + |y/|2 

Beltrami-Lapalace A dQ (f> = trV 2 <f> = V ■ V T </> = 1 ft (g ij y^l + \V f\ 2 dj( 

yi + lV/l 2 V 

where the matrix {g 11 ) = (8ij + difdjf)' 1 . For a C° n -fT 1 tensor T of (A;, Z)-type, the covariant 
derivatives T>T is a (fc, Z + l)-type tensor in L 2 . For any (k, I) tensor T\ and (k, I + 1) tensor T2 in 
C° n H 1 , one may verify, possibly through smooth approximations of d£l, 



[ (VT 1 )-T 2 dS= [ tr(Ta-r»T 2 (•,•)) ^. 

JdQ JdQ 



where, on the above right side, PT2(-, •) denotes the (fc,Z)-type tensor define by T>T2(X,Y)(- • • ) = 

(V X T 2 )(Y,y). 

From this identity, for any L 2 tensor T, one can define T>T, in the distribution sense, as in the 
dual space of C° D H 1 tensors. It is straightforward to verify that, for any (k, Z)-type tensor T in 
C° fl H 1 , we have 

(6.2) / T ■ A dn TdS = - [ \VT\ 2 dS. 

Jan Jdn 
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If T, VT G C° n H 1 , we have 



/ \A dn T\ 2 dS=[ \V 2 T\ 2 +\v n{x x )x T-K(X n ,X J2 )V Xn T 
Jan Jan 1 1 2 1 



• Vx T 



(6.3) - - \R,(X h ,Xj 2 )T\ 2 dS. 



i« 

Here {Xi, . . . , X n _i} is any pointwise orthonormal frame of 9fi which always appears in the trace 
form resulting in the independence of the corresponding quantities of the choice of the frame. The 
curvature acts on T in the usual sense 

R(X, Y)T = V Y V X T - V X V Y T - V [Y:X ]T = V 2 T(Y, X) - V 2 T(X, Y). 

Though TZ(X, Y)T seems to contain derivatives of T, one may calculate 

(K(X, Y)T)(X 1 ,X 2 , ...) = -T(K(X, Y)X 1 ,X 2 , ■ ■ ■) - T{X u Tl{X, Y)X 2 , ...)-.... 

So the dependence of 1Z(X, Y)T on X, Y, T is only pointwise and 1Z vanishes if T is a scalar function. 

Since / — Aqq is a positive self-adjoint operator on L 2 (dQ), for <f> : dfl — > K and r > 0, we define 
the Sobolev norm | • \u^(dU) on the surface dQ as 

\ 2 H r (dQ) = \(I ~ A dn )2(j)\ 2 dS; | • 1 1,2 (an) = I • \H°{dQ)- 

Jan 

As usual, for r > 0, | ■ \H~ r (dQ) coincides with | ■ |jj r (an)*- O ne may note here, since the Christofell 
symbols Tfj are only in L 2 , while \T\ h i^q^ < oo implies that VT £ L 2 (dU) from (|6.2[) . it does not 
imply that T is in H l in local coordinates, except when T is of (0, 0)-type, i.e. a scalar function. 
Similarly, from (fOjl. |T| H 2 (9n) < cxo does not imply £> 2 T € L 2 (dQ,). 

Remark. When n = 2, dVt is 1-dimension and Aqq = d ss where s is the arc length parameter, 
which is well defined if Oft is in W 1,1 . In fact, dQ G H s , s > \, is sufficient for the definitions of 
all the objects intrinsic in d£l. 

H s boundary (917, s > ^9— ■ For the purpose of this paper, we assume dtl C K" is in H s with 
s > ■ On the one hand, we defined the norm | • \H r (dfi) using the Beltrami-Lapalace Aan m the 
above. On the other hand, an obvious and traditional way to define the Sobolev space H r (d£l), 
— s<r<s for scalar valued functions and 1 — s < r < s — 1 for tensor valued functions, is through 
local coordinate coverings of dfl and the definition of the Sobolev space H r (R n_1 ). From standard 
Sobolev inequalities, it is easy to see that the latter definition of the spaces H r (d£l) is actually 
independent of local coordinates and naturally induces a topology on H r (dO,). In particular, when 
r > is an integer, straightforward calculation also shows that a function / (or tensor field) belongs 
to H r (dQ) if and only if f,V r f G L 2 (dn). In fact, we have 

Proposition 6.1. For r G [— s, s] (r G [1 — s, s — 1] for tensors), the norm \ ■ \n r (dQ) ^ s equivalent 
to the norm on H r (d£l) defined by using local coordinates. 

The proof of this proposition follows from the standard elliptic estimates using the local co- 
ordinates along with interpolation. In particular, when s is an integer, one may also prove it 
geometrically. In fact, the proposition clearly holds for r = 1 and r = 2 due to (|6.2|) and ()6.3jl and 
Sobolev inequalities. When r is an integer and r G [3, s] (r G [3, s — 1] for tensors), the proposition 
can be proved by using the following identity 

(6.4) ([A dn ,V]T)(X) = K(X,X j )Vx j T+<p x .K)(X,X j )T + (H(X,X j )V T){X 3 ). 

Finally, for non-integer or negative r, the proposition follows from interpolation and duality. An- 
other implication of (|6.4|) is that V : H r (dQ) — ► H r ~ 1 (dQ) is bounded and | • \L(H r (dn),H r - 1 (dQ)) 
depends only on 1Z and its derivatives. 

It is well known that the regularity of dfl can be determined from the regularity of its mean 
curvature k. 
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Proposition 6.2. Let Q C R n be a domain such that 8ft G H s ° , sq > Suppose |k| #3-2(90,) < 

oo wii/i s > so, i/ien d£l G -fP. 

Proposition 16.21 can be proved by using local coordinates and standard quasilinear estimates. 
Another proof can be based on the following identity which is also used in a priori estimates. 
Intuitively, let e : d£l — > W 1 be the imbedding, then II = —N ■ T> 2 e and k = —N ■ Aqq e, where II is 
viewed as a symmetric quadratic form. Thus it is expected that the difference AqqII — T> 2 k should 
be of lower order terms only. In fact, 

(6.5) -A an n = -p 2 K + |n| 2 n-Kn 2 . 

To prove the identity, at any x G let Aj, i = 1, . . . , n — 1, be the eigenvalues of n(x) and X{ 
be the associated eigenvectors which form an orthonormal frame of T x dVt. Parallel transport this 
frame to every base point in a neighborhood of x on dVt along the radial geodesies emitting from 
x. From the construction, we have T> x Xj = 0, [Xj,Jfj] = 0, and U(Xi) = A^Xj at x. For any 
X = a^Xj with constants a , . . . , a™ -1 , at x, 

(A 9n U)(X,X) ={V Xi V x Ji)(X,X) = V Xi ((V Xi IL)(X,X)) = V X ((V X U)(X U X)) 
={V Xi V x Ii)(X u X) = (V x V Xt U)(X,Xi) + (K&XiWiX^X) 
=V X ((V X .U)(X, Xi)) + U(TZ(X l ,X)X u X) + U(X i ,TZ(X i ,X)X). 

For the first term at x, from the construction of our special frame, 

VxUVxJtyfrXi)) = V x ((V x Il)(Xi,Xi)) = V x V x k - 2U(X i ,V x V x X i ) = V 2 k{X,X). 

To calculate the remaining two terms, one may substitute X = a J Xj and use H(Xj) = \jXj, the 
symmetry of 1Z, and the following calculation 

4JZ(Xi, Xj^Xi ■ Xj 2 =lZ(Xi, Xj 1 + Xj 2 )Xi ■ (Xj 1 + Xj 2 ) — lZ(Xi,Xj 1 — Xj 2 )Xi ■ (Xj 1 — Xj 2 ) 

=n(Xj, Xi)n(Xj 1 + Xj 2 ,Xj 1 + Xj 2 ) — u(Xi,Xj 1 + x,- 2 ) 2 
— n(Xj, Xi)n(Xj 1 — Xj 2 ,Xj 1 — Xj 2 ) + n(Xj, x^ — Xj 2 ) 2 

Equality (|6.5|) follows consequently. 

H s ° neighborhoods of domains, sq > Given a domain with dQ* in H s °, we will consider 
the set A = A(0*, Sq, s, L, 5), s > so, of neighboring domains of O*, given in Definition I4.2I From 
Proposition ^, 21 every £1 G A is in H s . Given and sufficiently small 5 > 0, in the following, we 
will derive some estimates with bounds C uniform in Q G A. Since dfi* is compact, for any a > 0, 
there exist xi G 1" and d, d% G (0, 5], i = 1, . . . , m, 

(Bl) B(dn*,d) C U^ 1 R l (d l ) where each i^(-) = x J 4 (.) C K n with ^(-) and !*(•) being an 
open (n — l)-dimensional disk and an open perpendicular segment in R n , both centered at 
Xi and of the given radius and half length, respectively; 

(B2) For each i, z = (z±, . . . , z n _\, z n ) = (z, z n ) being an Euclidean coordinate system on Ri(-) x 
Ii(-), there exists an H s function /*j : Ri{2di) — > Jj, so that 

(6.6) |/«b° < vdi, \ D f*i\c° < a i and ^* DRi(2di) = {z n > f*i(z)}. 

For any cr > with a fixed coordinate covering {Ri(di)}VL 1 of 9f2* of the above type, it is clear 
that, when S > is sufficiently small, {Ri{di)} 1 £ =1 is still a coordinate covering of any d£l G A 
satisfying (Bl) and (B2) with coordinate functions {/, G H s } r j ^ =1 . This will provide us some 
technical convenience in deriving estimates uniform in O, G A. The following proposition is a 
refinement of Proposition 16.21 
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Proposition 6.3. Given f2*, there exists 5 > such that, for any L > 0, there exists C > swc/i 
that the second fundamental form of any Q € A satisfies 

|n|H«-2(an) < C. 

Proof. The proof follows simply from the standard elliptic estimates and we will only give a sketch. 
With dfl £ H s due to Proposition \V>.2\ we will use the above coordinate covering {Ri(di)}™ =1 and 
the coordinate functions {fi G il s (i?j(2dj))}|^ 1; whose H s °{Ri(2di)) norms are uniformly bounded 
in i and S7. Let 7 : [0, +00) — > [0, 1] be a smooth cut-off function supported on [0, |] and j\^ Q 5] = 1. 

On each Ri(2di), let 



7i(^) = 7(^7)) = 7i(£) K (5, 9i = lift, 

where k is the mean curvature of 90. It is clear from the definition of A that \ K i\H s i- 2 (R-(d-)) * s 
bounded uniformly in i and O, for si = min{sg + 2, s}. From the mean curvature formula (|6.1|) . 



dj 1 fjdj 2 fj ^ _ n | [r7 ^|2^ 
1 + |V/,| 2 



-A 5i + ,^,; a ili2 ft =(1 + \Vfi\ 2 )2Ki - Ajifi - 2D 1% • Dfi 



@jl f$J2 fi 



( d jlhlifi + d hH d hfi + 9 hli d jji) 



i + fv/,1 2 

Since 7 is supported on [0, |], without loss of generality, we may treat fi as compactly supported on 
the ball of radius because fi can always be replaced by l(^p)fi(z) where 7 is a cut-off function 
supported on [0, |] and 7|[q 3] = 1- By partition of the unity and the Inverse Function Theorem, 
• • • > fm) can be expressed by g = (g±, . . . ,g m ) with the same regularity and similar estimates. 
Thus, dividing both sides of the above equation by (1 + |V/j| 2 )2 , it can be rewritten as 

-Ag + A jlj2 (x,g,dg)d jlj2 g = K + G{x,g,dg), 

where Aj 1 j 2 , G\, G2 are smooth in their arguments and Aj 1 j 2 < Co~ 2 so the left side is uniformly 
elliptic. In this form, the estimate on \g\H s i(M. n - 1 ) uniform in Q is obtained following the standard 
theory of quasilinear elliptic equations. If s± < s, this procedure can be carried out again with so 
replaced by si = sq + 2. Thus the desired uniform estimates on fi in H s follow by repeating this 
procedure. □ 

Remark. 1 ) One can also prove Proposition \6.tft based on (j6.5j) . 

2) A more careful estimate can be found in Lemma \4-2\ when a more detailed relationship between 
\H\fjs-2^QQ^ and \k\jjs-2(qq-) is given under certain conditions. 

Using (|6.2j) . (|6.3I) . (|6.4fl . and the above uniform estimate on II, which implies the uniform estimate 
on the curvature TZ, it is easy to prove that, for any tensor T € H r (d£l), r S [2 — s, s — 1] (r G [1 — s, s] 
for scalars), we have 

(6-7) \^T\H r -^(dn) < C\T\ H r( dn ) 

for some C uniform in f2 € A. 

From the uniform estimates on those (uniformly fixed) local coordinates derived in the above 
proof, it is also clear that the constants in the Sobolev inequalities (e.g. | • \H s (dn) to LP(d£l) or 
C a {dQ) for s < k) are uniform in O € A. The two most used inequalities in this paper are for 
/ G H Sl (dQ) and g G H S2 {dn), Sl < s 2 , 

Tl — 1 

\fg\ HS1+S2 -n^l „ n) < C\f\ H s l(dn) \g\ H s 2{dn) , if s 2 < — ^— and < si + s 2 

Tl — 1 

\fg\H*i(dn) < c \f\H s i(dn)\g\H s 2(dn), if s 2 > — ^— and < si + s 2 . 

Similar inequalities hold for / and g defined in S7. 
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6.2. Dirichlet-Neumann operator. Given fi*, in order to study the Dirichlet-Neumann opera- 
tor for domains fi£ A = A(fi*, sq, s, L, 5), we need to first construct local coordinate maps on each 
Ri(2di) for each fi, which flatten <9fi and have estimates uniform in € A, based on the above 
coordinates functions of dfi. 



Local coordinates and partition of the unit. From Proposition l6.3l d£lnRi(2di) is represented 
as the graph of an H s function fi : Ri(2di) — > R. Let (fi = 7j/j with ji defined in the previous 
proof. A standard way to extend (fi to a function € H s+ ^(M. n ) is through the Fourier transform 
with an appropriate constant a: 

(6.8) = a (1 + (gl)2 + --- + (r " 1)2) ^ (e 1 , • • • , f- 1 ) 

(i + (e i ) 2 + ... + (e) 2 r + ^ 

Since s + i > § + 1 and <1> is bounded in H s+ ? uniformly in i and fi, |-D<3?|<70 is also uniformly 
bounded. Therefore, there exists b > so that 

(6.9) H t (z\ ...,z n ) = {z 1 ,..., z n ~\bz n + $(z\ . . . , z n )) 

is a diffeomorphism so that \HA „, i and „, i are bounded uniformly in i and 0. 

Let Gj = (i?i) _1 and gi(z) be the n-the component of d, then |gi|^ s+ i (Rn) , d^ft, and (d^g^ 1 

are bounded uniformly in i and 0. Obviously, there exists a uniform 5* > so that 

5 5 5 5 

(Ri(-di) x Ji(-Mi)) n fi = /ii(-di) x /<(-5*di)) n {g, > 0}. 

Based on the local coordinate maps, we can construct partition of the unit satisfying estimates 
uniform in fi if 5 is small. In fact, take 7, £ : C°°([0, +00), [0, 1]) so that supp(7) C [0, |], j\^ Q 9j = 1, 

£(r) = r for r > |, and £|[o I] = 3- Define 

7„( 2)=7 (M) 7( M ); ^(oE^feoG,), 7 „ = 2^i, 7 * = (l-£™ l7 «Mfi). 
dj o*dj 77 

It is straight forward to verify that 7*0,7*1, • • • >7*m £ H s+ 2 (R™, [0, 1]) satisfy 

5 5 

bwljyH-J (R «) ^ C > supp(7*i) C flj(^di) x J^-^di)), 
for i = 1, . . . , m, and (E™ 7*;) \ n = 1. 



Remark. Using the above local coordinates and partition of unity we can establish the equivalence 
of the standard H e norm and the norm given in definition \6.1\ for integer £ G (— s — 5, s + 5)- ^ e 
raiio 0/ i/ie too norms is bounded above and below by two constants depending only on A. 

Trace and Harmonic extension. Let s\ £ (i,s + i]. Using the partition of the unit and the 
above local coordinates, it is straight forward to obtain the trace operator estimate 

(6-10) \man)\ Hn - km) <C\*\ H s lin) 

for any \& € H Sl (fi) where C > is uniform in fi € A. 

In order to obtain the estimate on the Harmonic extension operator, we first construct an ex- 
tension for convenience. Let S2 £ (0, s] and ifi € H S2 (d£l). Take the same auxiliary functions 7 and 
£ used above. For each 1 < i < m, let fa{z) = j(^-)tfi(Hi(z,0)) and $i(z) be the extension of fa, 
defined in the way of (|6.8j) . Let 

*i(z) = $ i ( Z ) 7 (|l) 7 (M ), * x = ^ o G< 

(Jjq (J ^ (Jl/£ 
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where can be viewed as a function defined on W 1 D £1- Let 

m(z) = 7(^) 2 7(^)> V = e o E^fa o Gi ), * = ^. 

It is easy to verify that ^ € H S2+ 2 (W 1 ) is an extension of ^ G H S2 (d£l) satisfying the estimate 

(6- 11 ) m H s 2+knn) <C\rP\ HS2m 

with C uniform in f2 € A. 

Using the partition of the unit and the local coordinates we constructed above and following the 
standard procedure, we have 

Lemma 6.4. There exists C > which depends only on the set A so that, for s\ G [o> s ] 

1^ — ll ^ I l^yl < C 

1 'l(H s i-1 (f2),/i" s i + 2(0)) 1 'L(_H" s l(f2),/i" s i + ^(Q)) _ 

Dirichlet-Neumann operator. Following from the above estimate, the Dirichlet-Neumann op- 
erator TV : H Sl (d£l) — > H Sl ~ 1 (dfl) can be defined and it has a uniform bound for s± S (1, s]. In 
fact, we can extend M into a weaker form defined on H Sl {d£l) for s\ > ^. Given / G H^{dVl), 
define 7V(/) € H-^(dQ) as 



<iP,Af(f) >= [ Vfn-V^ H dx 
Jn 



for any ^ € H 2 (017). It is easy to prove that 

1) M is self-adjoint in L 2 (dQ) with compact resolvent; 

2) the kernel ker(AA) = {const}; 

3) Cl/I , >\M(f)\ 1 > Al/I 1 for any / satisfying L n fdS = 0. 

The first inequality of 3) follows from the uniform bound on ft. In order to prove the second 
inequality in (3), one notices that 



From the estimate of the trace operator, we only need, for any / satisfying J dn fdS = 0, 

\n(f)\ L 2 {n) < c\vn(f)\ L 2 (n) 

with a constant C uniform in / and £1. This inequality can be proved by a compactness argument. 
Thus, by duality and interpolation, N can be extended to H Sl (dft) for all s € [1 — s,s] and 
\N\i,{H a i (dn),H s -i~ 1 (dn)) i s bounded uniformly in £1. Moreover, for / G H s (dQ) with Jg^fdS = 0, 
we can obtain — C|^(/)li? s - 1 (9n) with C uniform in fi. The proof is simply the elliptic 

estimate under the Neumann boundary condition - very much similar to the derivation of the 
harmonic extension estimate, except in the first step, instead of using (|6.8j) . we need to construct F 
with \F\ s+ i < C\ftf(f)\ HS -i(QQ\ and VatF = M(f) on Oft, by using a slightly different formula 

of the same fashion. Therefore, from interpolation, we have, for any s\ € [^,s], 

l/lfln (a!1 )<c|Ar(/)| ff ,-i (M)) if / fds = o 

with C uniform in f2. moreover, this inequality holds for s± G [1 — s, s] by duality. Based on these 
estimates, we can use / + AT to define the Sobolev norms which are equivalent to those defined by 
using / — Aqq uniformly in f2, i.e. 



30 SHATAH AND ZENG 

Proposition 6.5. For s\ G [— s,s], the norms on H Sl (dQ) defined by interpolating I — A T and 
I + J\f are equivalent, i.e. 

I(I-A af # < (I + M) S1 <C(I-A dn )% 

with C uniform in 0, £ A. 

Furthermore, for si G [— s, s — 1], 

AT 1 : H Sl (dn) -> H Sl (dn) = {/ G H Sl (dn) | / /dS = 0} 

Jan 

is well defined and bounded uniformly in $7. AA _1 denned on -H~2(dfl) induces the solvability of 
the Lapalace equation with Neumann boundary data given in H~2(dQ). 

To demonstrate that M behaves like differentiation, we give the following "product rule" . Given 
functions / and g defined on d£l. Since 

fn9n ~ A~ 1 A(f n g n ) = H(fngn\dn) = Uifg) in Q, 

we obtain 

(6.12) M(fg) = fM(g) + gAf(f) - 2V N A-\Vf H • Vgn)- 

Since M is like differentiation, coordinate independent, and self-adjoint, appearing naturally in the 
Euler's equation, it is sometimes convenient to express the Sobolev norms on dVt by A^. 

Relationship between J\f and Ago,. In addition to just the comparison between the norms 
of Aqq and ftf, we will prove that J\f is simply equal to (— Aqq)2 plus lower order terms. This 
improves the previous estimates and makes the estimates of some Sobolev norms using J\f more 
convenient. From the identity 

(6.13) Aip = A dn i; + KV N 4> + D 2 4>(N,N) x G dVl 

for any smooth function ip on VL. Recall that Nfi(x) and k-h(x), x € fi, denote the harmonic 
extension of the unit outward normal vector and the mean curvature of d£l. Given smooth / : 
dQ — » M, at any x G <9f2, 

D 2 f H (N,N) =V N V NH f H - Vf H ■ V N N n 

=V N (H((V N fn)\aa) + (-A)- 1 (-A)(V 7VW /w)) - M(N) ■ (M(f)N + V T /) 

=N\f) - 2V N (-A)-\DN n ■ D 2 f n ) - M(N) • (Af(f)N + V T /) 

which implies 

(6.14) (-Aan - M 2 )f = KAf(f) - 2V N (-A)-\DN n ■ D 2 f n ) - Af(N) ■ (Af(f)N + V T /). 
Proposition 6.6. 1) For s > ^^p, there exists C > uniform in O E A such that we have 

|A 9f2 + N 2 \ L ( H si (m)tH si-i( m)) <C, s' G [2 - s,s - 1]. 

2) For s G ' hIf) an< ^ s > 2, there exists C > uniform in G A such that we have 

n + 1 

\A dn +Af 2 \ , , „ +5+ < C, s'g(2-s, ). 

Proof. For s' > ^^y^ — s, the above inequalities follow directly from (|6.14j) and the estimates on Ti 
and N . Thus, by duality and interpolation, we only need to consider s' = \ or s' = ^ in each case, 
respectively. Let /, g : d£l — > R be smooth and harmonically extend into f2. Equality (|6.14|) yields 

f g(-A 9n - N 2 )fdS = [ ngAf (/) + gM(N) ■ (M(f)N + V T f)dS - 2 / DN H (Vg H ) ■ Vf n dx, 
Jan Jan Jn 
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which is sufficient to establish the estimate. □ 

Corollary. The proposition implies the commutator estimates 

(6.15) \[ A dn,^]\ L{{H s' [dnhH s'-2 {dn)) <C, J €[3-8,s-l] 



if s > and 



n + 1 

(6.16) |[Aan, A/11 , , <C, s'e(3-s, ), 

i/a G (^f 1 , 11 ^) ands> 2. 

We need the following abstract result for a more careful estimate on J\f. 

Proposition 6.7. Let X be a Hilbert space and A and B be (possibly unbounded) self-adjoint 
positive operators on X so that A~ 1 B and AB -1 are bounded. Suppose K = A 2 — B 2 satisfies that 
KB~ a is bounded with a G [0,2), then (A — B)B 1 ~ a is bounded. 

Proof. Let R = A - B. Calculating (B + R) 2 = B 2 + K, we obtain 

-BR - RB = R 2 - K, 

which implies 



Therefore, 



?L(e- Bt Re- Bt ) = e~ B \R 2 - K)e~ Bt > - e ' Bt K e - Bt > -Ce~ Bt B a e~ Bt . 

R<d e- Bt B a e- Bt dt = —B 01 ^ 1 
./n 2 



Calculating A = (A — R) + K with a similar procedure, we obtain 



Co 



R> -C 2 I e- At A a e~ M dt = -—A a ~ l . 
' Jo 2 

Thus, the conclusion follows. □ 

From Proposition 16.61 and Proposition 16.71 we obtain 
Theorem 6.1. There exist C > 0, which depends only on the set A such that if s > 

|(-A an )l -N\ L{H s> m)) <C, S '€ [1-8,5-1] 
and ifs>2ands€ (^±±, =±5), fora=^-s, 

1 Tl -\-~ 1 

(6.17) |(-Aai) 3 -N\ L{HS > im)tHs r- a+ i (m)) <C, S € (1-8,—^-). 

Proof. We will give the proof for the second case only as the proof for the first proof is similar. 
The estimate ()6.17j) follows directly from Proposition 16.61 and 16.71 for s' € (1 — s, ^jp). To prove 
the estimate for s' G [^r~, ^j^), we observe that s' — 2 G (1 — s, and we have 

|(I - A m y\(-A dn )l - *an)\ L{H ^ m) ^_ a+1{m)) < C. 

Thus (|6.17|) follows from the commutator estimate (|6.16j) . □ 

Decomposition of vector fields. We conclude this section by introducing the velocity field 
decomposition. Given an L 2 vector field u : 17 — > W 1 , it is standard to decompose it into the 
divergence free part v G L 2 and the gradient part —Vp for p G Hq(Q). In fact, 

(6.18) -Ap = V ■ u v = u + Vp. 
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For any divergence free vector field v G L 2 (Q), the normal component on the boundary v 1 - = v ■ N : 
dQ — > E in H~2(d£l) is defined as 

<t>~ L ,V>>= / v-V^-ftdx 
Jn 

for any ^ G Hz(dQ). By interpolation, for any si G [0, s — 5], 
(6-19) l^lfln-i^^^HH-itn), 

with C uniform in € A. This induces a decomposition of v into two divergence free parts, the 
rotation part v r and the irrotational (or gradient) part follows 

(6.20) t> ir = VHN^v- 1 -, Vr = v - Vir . 

It is easy to verify v r , Vi r G L 2 (Q) and 

V • v r = V • t>j r = 0, < v r , Vi r >= 0, v^r = 0. 

If v is a divergence free velocity field, v r component is responsible of the internal rotation and Vi r 
of the motion of the domain. 



Notation 



tr(.A): the trace of an operator. 

A*: the adjoint operator of an operator. 

Ai ■ A 2 = tr(Ai(A 2 )*), for two operators. 

B(S, e) = U X £sB(x, e): an e-neighborhood of a set S. 

D and d: differentiation with respect to spatial variables. 

V/: the gradient vector of a scalar function /. 

Vx- the directional directive in the direction X. 

_L and T: the normal and the tangential components of the relevant quantities. 

Dt = dt + v l d x i the material derivative along the particle path. 

D^: the projection of T) t to the tangent space of dVtt C E n . 

N(t, x): the outward unit normal vector of <9f2t at x G d£lf 

II: the second fundamental form of dflt, H(t,x)(w) = V W N G T x dQf 

U(X,Y) = U(X) Y. 

k: the mean curvature of <9$7t, i.e. K = trlT. 

jyi = H(f): the harmonic extension of / on f^. 

N(f) = Vat7^(/) : dft — > E: the Dirichlet-Neumann operator. 

X = X o u^ 1 the Lagrangian coordinates description of X. 

T>: the covariant differentiation on d£lt C E n . 

T^w = VjJ, for any x G dQ.t w G T x d£lf 

1Z(X,Y), 1,7 6 T x dCl t : the curvature tensor of dflf 

A_m = trZ? 2 : the Beltrami-Lapalace operator on a Riemannian manifold M. 

A -1 : the inverse Laplacian with zero Dirichlet data. 

r = {(f) : Q,t — > E n ; volume preserving homeomorphism} 

Q): the covariant derivative on T, 

<2): represent # in Eulerian coordinates. 

the curvature operator on T. 
M: represent in Eulerian coordinates. 
II: the second fundamental form of T C L 2 
ll u (w 1 ,w 2 ) = V^^, for any u G T, ^1,^2 G T„r 
Pv,w = — A _1 tr(Df Dw). 
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